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ABSTRACT 

We present Spitzer 3.6-24 /im photometry and spectroscopy for stars in the 1-3 Myr-old Coronet 
Cluster, expanding upon the survey of Sicilia-Aguilar et al. (2008). Using sophisticated radiative 
transfer models, we analyze these new data and those from Sicilia-Aguilar et al. (2008) to identify 
disks with evidence for substantial dust evolution consistent with disk clearing: transitional disks. We 
then analyze data in Taurus and others young clusters - IC 348, NGC 2362, and T) Cha - to constrain 
the transitional disk frequency as a function of time. Our analysis confirms previous results finding 
evidence for two types of transitional disks - those with inner holes and those that are homologously 
depleted. The percentage of disks in the transitional phase increases from ~ 15-20% at 1-2 Myr 
to > 50% at 5-8 Myr; the mean transitional disk lifetime is closer to ~ 1 Myr than 0.1-0.5 Myr, 
consistent with previous studies by Currie et al. (2009) and Sicilia-Aguilar et al. (2009). In the 
Coronet Cluster and IC 348, transitional disks are more numerous for very low-mass M3-M6 stars 
than for more massive K5-M2 stars, while Taurus lacks a strong spectral type- dependent frequency. 
Assuming standard values for the gas-to-dust ratio and other disk properties, the lower limit for the 
masses of optically-thick primordial disks is M<n s k ~ 0.001-0.003 M*. We find that single color-color 
diagrams do not by themselves uniquely identify transitional disks or primordial disks. Full SED 
modeling is required to accurately assess disk evolution for individual sources and inform statistical 
estimates of the transitional disk population in large samples using mid-IR colors. 

Subject headings: 



1. INTRODUCTION 

Transitional protoplanetary disks identify important 
stages in planet formation, linking epochs when disks 
are optically-thick primordial protoplanetary disks, rich 
in gas and containing massive amounts of small dust, 
and when they are gas-poor/free, optically-thin debris 
disks and thus are in the post gas pl anet-building stage 
(|Currie et al.ll2009HStrom et al.lll989l) . Compared to pri- 
mordial disks, transitional disks have weaker near-to-mid 
IR emission indicating that they are depleted in warm, 
micron-sized dust within several AU of the host star (e.g. 
iStrom et al.|[l989HCalvet et al.|[200l . If this dust traces 
the bulk mass of solids, transitional disks then show ev- 
idence for the removal of planet-building material. 

Though transitional disks are not primarily identi- 
fied by their levels of gas, they ty pically have a lower 
rate/frequency of stellar accretion (ISicilia-Aguilar et al.l 
I2006at iMvTzerolle et al.ll20H iCieza et al.ll20100 . Transi- 
tional disks may also be losing much of their gas from 
a ph otoevaporative wind driven by stellar EUV pho- 
tons dPascucci and Ster zik 2009), may lack gas-rich in- 
ner disks (Sal vk et al.ll2009f). or may have a low density 
of cold outer disk gas ([Hughes et al.|[2010l ). Thus, tran- 
sitional disks are also plausibly depleted of circumstellar 
gas responsible for the envelopes of gas/ice giant planets 
and may be signposts for the last phases of protoplane- 
tary disk evolution and gas giant planet formation. 

Spitzer studies have better constrained the morpholo- 
gies and lifetimes of transitional disks around solar to 
slightly subsolar-mass stars (M* ~ 0.5-1.4 M , see 
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ICurrid [20101 and [Muzcrollc 2010 for reviews). Transi- 
tional disks exhibit one of two general morphologies: 1) 
disks with weak or negligible levels of near-to-mid IR 
emission indicative of inner dust holes but more strongly 
emitting, typically optically-thick outer disks consistent 
with an inside-out clearing of dust or 2) disks lacking 
evidence for inner holes but having weaker emission at 
all IR wavelength s consistent with a more homologous 
depletion of dust (jCurrie et al. 2009). Transitional disks 
comprise < 1 0-20% of the disk popu l ation in « 1 Myr-old 
clust ers (e.g. iMuzerolle et al.l 120101 : ICurrie and Kenvonl 
2009) but are relatively more frequent by ~ 5 Myr, 
consistent with a transitional disk phase comprising a 
substantial fraction of the total proto planetary disk life- 
time (e.g. ~ 1 Myr out o f 3-5 Myr. (Currie et al.l 120091 
ISicilia-Aguilar et al.l 120091) 



Th e 1-3 Myr-old Coronet Cluster (|Mever and Wilkingj 
l2009l and references therein) may provide a particu- 
larly useful probe of transitional disk properties in the 
youngest clusters a n d as a function of stellar mass. 
ISicilia-Aguilar et al.l ([2008) find that transitional disks 
comprise nearly half of the disk population, much larger 
than frequen cies derived for other clus t ers of compara- 
ble a ge (e.g. ICurrie and Kenvonl 120091 : IMuzerolle et al.l 
l2010h . A t first glance, these r e sults appear to be in con- 
flict, but ISicilia-Aguilar et al.l ([2008T ) note that the mem- 
bers in_thejr_sj : m£k2jiave^ystematically low masses (see 
also iLopez-Marti et al.|[2010t ). Compared to solar and 
slightly subsolar-mass stars, disks around very low-mass 
stars and brown dwarfs may have very diff erent struc- 
tures and different dispersa l timescales (e.g. lLada et "all 
l2006t Hartman riet al.l l2006). Thus, it is unclear whether 
the high reported transitional disk frequency in the Coro- 
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net Cluster is anomalous or identifies a stellar mass- 
dependent transitional disk frequency. Addressing this 
issue requires analyzing disks around higher-mass Coro- 
net stars and assessing the transitional disk frequency as 
a function of stellar mass for many clusters. 

In this paper, we investigate the transitional disk pop- 
ulation in the Coronet Cluster as a function of spec- 
tral type/stellar mass by presenting Spitzer photome- 
try and spectroscopy for additional, predominantly solar- 
mass members. By comparing results for solar mass 
and subsolar-mass members, we determine whether the 
Coronet Cluster disk population shows evidence for a 
stellar mass-dependent frequency of transitional disks. 
§2 describes our new observations, data reduction, and 
photometry/spectroscopy. In §3, we model the SEDs of 
disk-bearing stars to determine plausible di sk evolution- 
ary s tates following previous methods (e.g. ICurrie et al.l 
2009). Finally, in §4 we compare these results to those 
for other clusters to investigate the duration of the tran- 
sitional disk phase. 



2. DATA 



2.1. 



Spitzer IRAC and MIPS Photometry 

To t he sample of stars analyzed by lSicilia-Aguilar et al.1 
lj2008f ). we add IRAC and MIPS photometry of other 
Coronet Cluster members located on the processed IRAC 
and MIPS mosaics obtained from Program IDs 6 (AORs 
3650816 and 3664640), 248 (AOR 13469696), and 30784 
(AOR 17672960). Image processing and photometry for 
th ese stars was performed fo llowing methods described 
in lSicilia-Aguilar et al.l (|2008| ). Briefly, the basic calibra- 
tion data (BCD) frames were mosaiced together using 
MOPEX using standard input parameters for the over- 
lap correction, pixel interpolation, and outlier rejection. 
Aperture photometry was performed with APEX using 
a 3 pixel aperture, a 12-20 pixel background radius, and 
aperture corrections listed in the IRAC data handbook 
(1.112, 1.113, 1.125, and 1.218 for the four channels, re- 
spectively). For MIPS, we used a 5 pixel aperture with 
a 8-13 pixel background annulus, and the aperture cor- 
rection listed in the MIPS data handbook (1.167). 

We select cluster s tars from the catalog of 
iForbrich and Preibischl (|2007f ). which identifies x- 
ray bright sources with near-IR/optical counterparts as 
Coro net Cluster members. We add 32 cluster members 
from IForbrich and P rcibisc h) (|2007f) to the list of stars 
analyzed in lSicilia-Aguilar et al.l (|2008l) . The IRAC and 
MIPS coverage areas do not completely overlap, and 
many cluster members lack photometry in two of the 
four IRAC channels ([3.6] and [5.8] or [4.5] and [8]) or 
lack MIPS data. Other sources, particularly Class I 
protostars and bright, higher-mass stars (e.g. R CrA) 
saturate the IRAC and/or MIPS detector. Table [T] lists 
the names and coordinates for these targets. Table [2] 
lists their photometry. Sources with an "FP" prefix 
for their names were classifi ed as members only by 
IForbrich and Preibischl (|2007f ): the numbers following 
this prefix correspond t o the row number in Table 2 of 
IForbrich and Preibischl (|2007| ). 

To further assess photometric quality, we inspected the 
processsed IRAC and MIPS mosaics and compared the 
known target positions to the computed target centroids. 
Most sources lie in regions of low background and are 
thus uncontaminated by nebulousity. Most sources lo- 



cated in high-background regions are much brighter than 
the background, which is typically uniform. However, 
two sources - TY CrA and HD 176386B - are located 
in regions of high and highly-variable background. In 
the longer wavelength 5.8 /im and 8 fim channels, the 
computed centroid positions are well offset from the ap- 
parent star positions, especially for HD 176386B. Figure 
[T] illustrates this offset for the 5.8 yjm channel. 

To investigate how the mid-IR fluxes from TY CrA and 
HD 176386 are affected by nebular emission, we recom- 
puted the flux for these sources using much smaller aper- 
tures and background annuli, ranging from our default 
assumptions to a small 2 pixel aperture radius and 2-6 
pixel background annulus. Depending on the choice of 
aperture radii and background annuli, their fluxes at 5.8 
\xm and 8 [im vary by ~ 50%-200%. Therefore, we con- 
sider the 5.8 lira and 8 [im photometry for these sources 
to be unreliable. Since the nebular background emission 
increases with wavelength, it also renders fluxes for these 
sources uncertain at 24-100 [im. 

2.2. Spitzer IRS Spectroscopy 

To supplement our photometric data, we add spectro- 
scopi c data from the Infrared Spectrograph (|Houck et al.1 
l2004h spectroscopy for several intermediate-mass mem- 
bers: S CrA, V709, T CrA, IRS7w, and IRS5. The spec- 
tra were processed starting from the BCD data and us- 
ing the spectral extraction tools developed for th e FEP S 
Spitzer science legacy team (jBouwman et al.1 [2 0081. 



based on the SMART software package (jHigdon et al.1 
2004). The spectra were extracted using 6 and 5 pixel 
wide apertures in the spatial dimension for short- (7.5 
- 14 /im) and the long-wavelength (14 - 35 fim) mod- 
ules, respectively. The background was subtracted using 
associated pairs of imaged spectra from the two nod- 
ded positions along the slit, also eliminating stray light 
contamination and anomalous dark currents. The long- 
wavelength modules were not available for S CrA, and 
in the cases of IRS7w, IRS5, and T CrA, strong back- 
ground contamination and/or extended emission makes 
the spectral extraction uncertain at A > 14/im. There- 
fore, we only include in the analysis the short-wavelength 
region for these sources. 

2.3. Ancillary Data 

To compare the stellar and circumstellar prop- 
erties of all Coronet Cluster members, we add 
IRAC /MIP S photometry and I RS sp ectroscopy pre- 
sented in ISicilia-Aguilar et "all ( 20081) ; optical pho- 



tomet ry from lLopez-Marti et aL 



(2005), iTorres et al.1 



(2006) and from the SIMBAD Astronomical Database 

(various source s); near-IR photometry from 2MASS 

ISkrutskie et~aT1 (|2006T) and iMever and Wilkind (|2009ft . 
Spectral types for most cluster stars wit h new 
Spi tzer data d eri ve fr om IMever and Wilkingl (2009) 
or iNisini et alj (120051). Most stars studied in 
ISicilia-Aguilar et al.l (|2008l ). have spectral types derived 
from VLT/FLAMES. 

We also consider new, higher signal-to-noise 
VLT/FLAMES spectra of Coronet Cluster mem- 
bers, which yield new spectral types for several low-mass 
cluster stars - G-l, CrA-4107, CrA-4109, CrA-432, 
CrA-468, CrA-452, G-85, and G-87. These data will be 
discussed and analyzed in detail in a later contribution 
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(Sicilia-Aguilar et al., 2011, in preparation). Briefly, the 
medium resolution (R ~ 6000) spectra were centered 
on three bands covering 0.6-0.9 /im (A c = 0.682, 0.773, 
and 0.881 /im). For spectral classification, we first 
determined whether the stars were likely earlier or later 
than M3 using standard spectral indices. For later stars, 
we derived spect ral types using spect r al indi ces of TiO 
bands following ISicilia-Aguilar et ail ((2008). Earlier 
stars were classified by c omparing their spect r a to C ep 
OB2 spectra described in ISicilia- Aguilar et al.l (|2005l) . 

We estimate the spectral types of stars without 
spectra by deriving the best-fit effective temperature 
from model th e opti cal to near-IR SED using the 
iRobitaille et al.l (l2006[) radia tive transfer grid and then 
using the lCurrie et al.l (|2010f) effective temperat ure scale 
for stars M2 or earlier and lLuhman et aTl (|2003l ) for later 
stars. For some sources, we cannot derive a spectral type 
because the source lacks optical data or is a protostar. 
Tables [3] and 0] list the spectral types for each star. 

To probe the mid-IR fluxes of bright stars that sat- 
urate the IRAC and MIPS detectors, we combined our 
source list with IRAS da t a for RCrA, TY CrA, and S 
CrA from iWilking et all ([i"985L Il992h . Our IRAC and 
MIPS mosaics indicate that the emission from R CrA 
and S CrA clearly dominates over the background nebu- 
lar emission and thus are only weakly contaminated. As 
mentioned previously, TY CrA is heavily contaminated 
by background emission. Given the large beam size of 
IRAS compared to Spitzer, the IRAS flux quoted for TY 
CrA is also likely unreliable. 

2.4. IRAC and MIPS Colors and Observed SEDs of 
Coronet Cluster Stars 

To provide a first-order investigation of the Coronet 
Cluster disk population, we analyze the distribution of 
near-to-mid IR colors and optical/IR spectral energy 
distributions (SEDs). The left panel of Figure [5] dis- 
plays the observed IRAC colors of cluster stars with 
new photometry (b l ack d ots) and stars analyzed by 
ISicilia-Aguilar et ail f|2008l grey dots). The distributions 
appear similar. Many stars concentrate near zero color 
(as expected for bare stellar photospheres) and around 
[3.6] -[4.5] = 0.3-0.6 and [5.8]-[8] = 0.5-0.7. Compared 
to the typical IRAC colors of primordial disks i n 1-2 
Myr-old clusters (solid box. lHartmann et al.ll2005l ). col- 
ors for Coronet Cluster stars may be slightly bluer by ~ 
0.2 mag. 

The longer wavelength colors for cluster stars exhibit 
a wider diversity (F igured right panel). At le ast three 
sources, all from the ISicilia- Aguilar et al.l (|2008f) sample, 
show evidence for weak/negligible 5.8 /im excess but a 
> 4 magnitude excess at 24 /im characteristic of transi- 
tional disks with inner holes. Others have 24 \xm and/or 
5.8 /im excess emission characteristic of more homolo- 
gously depleted transitional disks, su ch as those found i n 
NGC 2362 (dotted enclosed region. ICurrie et al.|[2009h . 
which is weaker than e mission for the lowe r-quartile me- 
dian Taurus SED from lFurlan et al.l (|2006l lower-left tri- 
angle) . The cluster also includes many stars with strong 
IRAC and MIPS excesses comparable to the median Tau- 
rus and upper-quartile Taurus SED. Because only six 
stars in the new sample have K s , [5.8] and [24] detec- 
tions, it is less clear from the K s -[5.8]/K s -[24] diagram 
whether the distribution of colors from the new cluster 



sample displays the same diversity. 

Examining the SEDs of new sources clearly shows that 
they have a wide range of morphologies (Figure[3]). Many 
sources (e.g. IRS2, IRS5ab, IRS7w) have rising near- 
to-mid IR SEDs indicative of Class I protostars. Well- 
known sources with optically-thick disks, such as the 
intermediate-mass stars S CrA and R CrA, have nearly 
flat SEDs from 1 /im to 8 /im. 

IRS spectra strengthens our identification of sources 
whose photometric data alone make them difficult to 
classify. In particular, V709 shows evidence for a weak 
excess based on its MIPS-24 /im flux and IRS spectra. 
Its SED gradually peels away from the stellar photo- 
sphere from ~ 8 /im to 30-35 /im. Based on their 10 
/im absorption features, IRS-5 an d IRS-7w are likely 
protostars with cool envelop es (e.g. iWillner et all 119821 ; 
lAndre and Montmerld I1994T). consis t ent with previous 
classifications dHenning etail fT99l iChini et al.1 l200l 
iGroppi et al.ll2007t ). 

3. ANALYSIS 

To assess the evolutionary states of Coronet Cluster 
dis ks, we analyze source SEDs from the new sample and 
the ISicilia-Aguilar et "aTl (j2008f ) sample using both sim- 
ple theoretical comparisons and sophisticated radiative 
transfer modeling. We follow a sl ightly modifi e d vers ion 
of the analysis methods used in ICurrie et al.l (j2009T ) to 
distinguish disks in different states. These classifications 
are based solely on 1-24 /im data: since far-IR/submm 
data provides crucial constraints on disk properties, in 
particular the inferred disk mass, we consider these clas- 
sifications to be provisional. We will revisit these classi- 
fications later after analyzing stars in other clusters with 
longer wavelength data (§4). 

In our classification, we consider protostars and disks 
with three main evolutionary states: primordial disks, 
transitional disks, and debris disk candidates. We iden- 
tify protostars as sources with rising SEDs from the near- 
IR to mid-IR consistent with emission predominantly 
from a cold dusty envelope (T < 200-300 K). Primor- 
dial disks are disks that are optically thick (tjr >> 1) 
and thus lack evidence for dust evolution consistent with 
active disk clearing. 

Transitional disks are disks that show substantial dust 
evolution consistent with active disk clearing. They can 
have one of two morphologies. First, they can exhibit 
a drop in the optical depth (r ~ 1 or less) of emitting 
dust either at all IR wavelengths consistent with a re- 
duced disk mass (a homologously depleted transitional 
disk). Second, they can also have photospheric or weak, 
optically-thin emission at shorter wavelengths (e.g. 3.6- 
5.8 /im) and more optically thick emission at longer wave- 
lengths consistent with inside-out disk clearing (a tran- 
sitional disk with an inner hole). In this l atter group, we 
inclu de so-called 'pre-transitional disks' (Espaillat et al. 
l2007f h These disks have a small amount of high opti- 
cal depth material very close to the star but otherwise 
have large, dust-poor inner regions (e.g. large gaps) and 
optically-thick outer regions and thus also identify an 
inside-out disk clearing. 

Debris disk candidates lack evidence for accretion 
and have very optically-thin emis sion {tjr << 1) 
from a low mass (< 0.01 M,„r, n r., IChen et al.1 12005) 
of second-generation dust (|Backman and Parescd 119931 : 
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ICurrie et al1l2008t iKenvon and Bromlevll2008D . Many of 
the Coronet Cluster stars are M dwarfs with very low 
masses. Even if these stars have tenuous emission and 
lack evidence for accretion, they also have grain removal 
timescales comparable to the host st ar's age. Thus, their 
dust need not be second generation (jCurrie and Kenvonl 
2009), they could be very late-stage transitional disks, 
and at best can only be considered debris disk candi- 
dates. Our basis for identifying sources with tenuous 
24 iim excess emission as debris disk candidates follows 
ICarpenter et al.l ()2009|) and is purely empirical. They 
identify debris disk candidates around late-type stars as 
those with 0.25-0.5 24 fim excesses. In our analysis, 
we consider the evolutionary state of these tenuous disks 
around Coronet Cluster stars as indeterminable and sim- 
ply label them as transitional disks/debris disk candi- 
dates 3 . 

For most sources, only photomet ric data is available. 
For s ources with IRS data from ISicilia-Aguilar et all 
we add the IRS flux densities, which when com- 
bined with our photometry yields flux densities at equally 
spaced wavelengths (Table [5]). The inclusion of the IRS 
data provides better sampling of the SED from 5 to 24 
/zm and extends the wavelength coverage. 

3.1. Method for Disk Identification 
3.1.1. Fiducial SEDs 

We first compare source SEDs to geometrically flat, 
optically-thick models to identify disks with evidence 
for a reduced optical depth of emitting material (tir 
< 1). Figure |4] shows two flattened reprocessing disk 
models appropriate for a T e = 3850 K primary (M0 
spectral type), though SEDs for earlier and later stars 
(e.g. K5 to M6) yield similar mid-IR fluxes relative 
to the stellar photosphere. The first model (dash-three 
dots/diamonds) is the s tandard, razor-thin flat repro- 
cessing disk model from | Kenvon and Hartmannl (1987, 
see also lAdams et al.l Il987| ). truncating the inner disk 
at a dust sublimation temperature of T = 1500K. 
The second model (thick grey dashed line) is produced 
from the Whitney Mo nte Carlo radiative trans fer code 
(jWhitnev et al.1 l2003aH b1: iRobitaille efall 12006ft . The 
model assumes a disk mass of 0.05 M Q , no flaring (H/r 
= constant), no accretion, and no protostellar envelope 
emission. The disk is optically-thick to its own radiation 
over spatial scales relevant for our study (< 10 AU). 

The Whitney disk mo del consistently has wea ker disk 
emission than the simple IKenvon and Hartmannl flat disk 
model, because stellar photons in the Whitney model are 
partially attenuated before they reach the r = 1 surface, 
due to the disk's finite thickness. Thus, the disk is colder 
at a given location and reradiates less energy. To provide 
an empirical comparison for source SEDs, we overplot the 
lower-quartil e Taur us SED (Figure HI left panel) from 
iFurlan et al.l ((2006), which has been used in other work 

3 While some debris disks have much stronger emission, nearly 
all of them surrou nd what are (or wi l l be) main sequence A stars or 
early F stars (e.g. Rickc et al. 2005; Curric et al. 2008) that have 
masses > 1.5 Mq. Debris disks detected around lower-mass stars 
(the focus of this paper) are far more infrequent and typi cally have 
24 )im excesses of only a few tenths of magnitudes (e.g. IPlavchanl 
120091) . Our results are insensitive as to whether or not stars with 
tenuous excesses are considered to have debris disks or transitional 
disks. 



to help distinguish disks in different states (ICieza et al.l 
120081 : ICurrie et al.l 120091 : ICurrie and Kenvonll2009D . The 
lower-quartile Taurus S ED almost perfectly tracks the 
IKenvon and Hartmannl flat disk model. 

In the right panel of Figure |4j we compare our 
fiducial SEDs to a flared optically-thick disk model 
(thick grey long-dashed li ne) with H/r = r 2 / 7 (e.g. 
iChiang and Goldreichlll997l) using the Whitney code and 
a disk that has undergone signific ant dust settling (t hick 
black long-dashed line). Following lLada et al.l (|2006| ). we 
quantify dust settling by reducing the disk scale height 
by a factor of 3 compared to the thermal equilibrium 
value adopted in the flared disk model 4 . The Whitney 
settled disk model predicts emission roughly equal to the 
Whitney flat disk model through ~ 7 fim and greater 
emission at longer wavelengths in spite of being flared: 
this probably happens because the puffed up inner wall 
effectively attenuates flux from disk regions slightly ex- 
terior (e.g. the disk only comes out of shadow at regions 
that effectively emit at A > 7 fim). 

In Figure we plot the K s -[5.8,8] vs. RV[24] colors 
for the flattened and settled disk models as a function of 
disk inclination, produced by convolving the model out- 
put with 2MASS and Spitzcr filter functions provided as 
a part of the Whitney code. There is little variation in 
color, especially for the flat model, except when the disk 
is viewed nearly edge on. In this case, the disk heav- 
ily extincts the star, causing its observed colors to be 
far redder than in face-on cases. The loci of colors indi- 
cate that flattened optically-thick disks should typically 
have K s -[5.8,8, 24] > 1.1, 1.5, 3.5. Settled disks have 
slightly bluer K s -[5.8] colors (~ 0.85-1) but redder K s - 
[8,24] colors (1.6-1.9). For the models where the disk is 
not viewed edge on, the disks have K-band excess emis- 
sion of ~ 0.1-0.2 magnitudes, consistent with a visual 
inspection of Figure |U Thus, optically-thick, flattened 
disks lacking K band excess (e.g. because of submicron- 
sized grain growth), should then typically have K s -[5.8] 
> 1.2 and K s -[8] > 1.6-1.7. The colors for an optically- 
thick flared disk lie off the plot range (e.g. K s -[8] = 
2.75, K s - [24] = 6.75). The disk emission at 24 fim varies 
wildly with disk flaring, so MIPS-24 fim does not effec- 
tively probe the disk optical depth. The colors from disk 
models for K5-M6 stars agree to within ~ 0.05-0.1 mag 

Since the Whitney flat disk model produces weaker 
emission through 8 fim than the other fiducial models, 
it defines a conservative limit for disks with optically- 
thick emission. Because we identify primordial disks as 
those with optically-thick IR emission, disks with 2-8 
fim emission greater than this flat disk model are consis- 
tent with being primordial disks. Conversely, disks with 
weaker emission correspond to later evolutionary stages: 
transitional disks and debris disk candidates. 

We emphasize that our adopted optically-thick disk 

4 Technically, we assume less "se ttling" than in other " settled 
disk" models, like those presented in[D'Alcssio ct al. (2006). How- 
ever, their treatment of settling assumes that some small, strongly 
emitting grains remain at larger heights above the disk midplane, 
whereas the Whitney code lacks this superheated grain population 
and thus may und erpredict near-to - mid I R disk emission. Further- 
more, some of the ID'A lcssio ct al. (2006) "settling" models likely 
produce weak emission because they assume low disk masses, in 
addition to being "settled". 
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limit is conservative. Emission in excess of this limit at 
A ~ 5-8 fim may be explained by disk flaring. While 
the model used to produce this limit lacks flaring, real 
disks almost assuredly have some flaring. Furthermore, 
the IRAC 8 /im channel overlaps with the fO p,m sili- 
cate featu re, complicating comp arisons with the flat disk 
limit (sec Muzcroll e et al J 1 201 0). The model also lacks 
accretion luminosity, which contributes some IR flux, but 
many disks we model show unambiguous accretion sig- 
natures. 

3.1.2. SED Modeling 

After comparing source SEDs to our flat repro- 
cessing disk SED, we model optical-to-mid IR fluxes 
with the grid of radi ative transfer disk models from 
iRobitaille et al.l (|2006|) . which yields estimates for disk 
properties (e.g. mass, i nner radius) consistent with other 
independent work (see Robitail le et al.l l2007). By iden- 
tifying the best- fitting models (min(x 2 )), we determine 
whether the disk plausibly has an inner region cleared 
of dust and is thus consistent with being a transitional 
disk with an inner hole. The computed inner disk radius 
distinguishes disks that lack inner holes (primordial, ho- 
mologously depleted transitional disk) and transitional 
disk with inner holes. To be conservative, we set the 
division at ten times the dust sublimation radius: disks 
with Ki n > 10 R su b are identified as transitional disks 
with inner holes. To account for flux calibration uncer- 
tainties, photometric errors, and variability, we assume a 
minimum 10% flux uncertainty in each photometric filter 
and IRS monochromatic flux density. 

We restrict ourselves to disk models which provide a 
good fit, which wc define as x 2 -X 2 best < 3, where x 2 best 
corresponds to the minimum \ 2 P er datapoint. This cri- 
teria is very similar to that used in previous studies (e.g. 
Robitaille et~aTll2007l : ICurrie et all 120091 : lErcolano etall 
2009). Finally, we incorporate systematic uncertainties 
in the extinction Ay and distance. For the extinction, 
we nominally assume a 20% systematic uncertainty with 
a "floor" of ± 0.5 mags. For the distance, we assume 150 
pc with an uncertainty of ± 20 pc. 

3.2. Results 

Tables [3] and [4] describe our provisional disk modeling 
results. The combined sample from both tables consists 
of 5 protostars, 16 stars with primordial disks, 6 with ho- 
mologously depleted transitional disks, 4 that have tran- 
sitional disks with inner holes, and 6 bare stellar photo- 
spheres. Sources not analyzed here lack sufficient photo- 
metric data to model their SEDs and/or determine their 
plausible range of spectral types (e.g. FP-8, FP-33, FP- 
37, and FP-38), making it impossible to identify their 
disk evolutionary states. Below we discuss modeling of 
several sources illustrating how we determine disk states. 

3.2.1. Disks in Different Evolutionary States 

Figure El co mpares the SED for CrA -159 from the new 
sample to the iKenvon and Hartmannl razor-thin flat disk 
model and the flat disk model produced from the Whit- 
ney radiative transfer code. The observed SED is dered- 
dened from Ay = 3. CrA-159 has an infrared excess de- 
tectable from 1.6-2 /jm to 24 /jm, which consistently lies 
well above the Whitney flat disk model, indicating that 



its emission is optically thick. According to our criteria, 
this source then has a primordial disk. With the ex- 
ception of one star (V709), all disks in the new Coronet 
sample, predominately comprised of higher-mass stars, 
yield emission lying above the flat, optically-thick disk 
limit. 

The disk population for the lSicilia-Aguilar et al.l (|2008l ) 
sample is far more diverse, including many transitional 
disks. We show comparisons between three homologously 
depleted transitional disk SEDs and flat disk models in 
Figure [7] Each source has a reduced level of disk emis- 
sion compared flat disk models. Likewise, all transitional 
disks with inner holes have weak/negligible 3.6-5.8 \xm 
emission but optically-thick 24 /im emission. 

The top panel of Figure[8]shows model fits for CrA-205, 
which has a transitional disk with an inner hole. CrA-205 
exhibits no excess emission shortwards of ~ 10 fim but 
has a rising SED clearly departing from the stellar pho- 
tosphere by ~ 15 [im. All of the best-fitting Robitaille 
models have inner disk radii (Rj„ > 1000 R s «b) much 
larger than our threshold identifying transitional disks 
with inner holes (10 R su &). Moreover, the grid sampling 
of inner disk radii (grey shaded region) is heavily peaked 
at 1 K su b (e.g. no inner hole). In spite of this intrin- 
sic bias against selecting disk models with holes, all disk 
models lacking inner holes fail our y 2 thre shold. 

Our results for the ISicilia-Aguilar et al.1 sample are in - 
ter mediate between those of ISicilia-Aguilar et all (2008) 
and lErcolano et al.1 (|2009l ) who disagreed over the frac- 
tion of transitional disks with inner holes. CrA-466 
and G-65 have best-fit models that do not require 
an i nner hole (Figure [H middle and bottom panels). 
Like lErcolano et al.l (|2009l ) , we identify them as primor- 
dia l disks. However, o u r modeling supports the claim 
by ISicilia-Aguilar et all (|2008l I2009D that CrA-205 has 
an inner hole, in fac t the largest one in our sample. 
lErcolano et al.1 (|2009h did not to identify CrA-205 's in- 
ner hole probably because they do not include the IRAC 
8 fim data nor any IRS dat a in their fitting . 

While we also agree with lErcolano et al.l ([2009) that 
G-14 and G-87 lack evidence for an inner hole, our mod- 
eling suggests that both have weaker emission than a per- 
fectly flat, optically- t hick reprocessing disk. Thus, like 
ISicilia-Aguilar et "all (|2008| ) we identify these as transi- 
tion al disks, albeit ones th at are homologously depleted. 
The lErcolano et ah! (|2009l ) study focused on identifying 
transitional disks with inner holes. Therefore, we identify 
more transitional disks because we adopt a more expan- 
sive definition for what constitutes a transitional disk. 

Two of the sources labeled as homologously depleted 
transitional disks may be debris disk candidates, while 
one of the transitional disks with inner holes may be a de- 
bris disk candidate. The identity of V709 is particularly 
questionable: while we identify it as a disk-bearing star, 
it is possible that its weak excess could be due to nebular 
contamination or an offset between the near-IR data and 
the MIPS data/IRS spectra caused by variability such 
as that produced by an undetected eclipsing companion. 
More importantly, if it has a disk, it is unclear whether 
it is a debris disk or transitional disk. 

Cluster members with new Spitzer photometry are 
mostly protostars or stars with primordial disks (5 and 8, 
respectively, Table [H). The disk popula t ion fo r M stars 
sample studied in Sici lia-Aguilar et all (|2008l ) is more 
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evenly divided between primordial disks (8) and tran- 
sitional disks (3-4 with inner holes, 4-6 homologously 
depleted). Unlike t h e new sample, many stars in the 
ISicilia-Aguilar et "all ((2008) sample have homologously 
depleted transitional disks. Thus, based on modeling 
1-24/70 fim data, the percentage of protoplanetary (pri- 
mordial + transitional) disks in our combined sample 
that are transitional disks is ~ 30% (7/23). The percent- 
age rises to 41% if debris disk candidates are classified as 
transitional disks, including V709 (11/27). 

3.2.2. Disk Properties As a Function of Stellar Mass 

To examine the spectral type/stellar mass dependence 
of disk properties, we bin our sample into three spectral 
type groups: earli er than K5, K5-M2 , and later than 
M2. Assuming the Baraffc c t al.1 (|1998| ) isochrones (mix- 
ing l ength = 1.9 H P ) an d the T e vs. spectral type scale 
from lCurrie et all (|2010D . the corresponding stellar mass 
division is M* > 1 M , = 0.5-1 M Q , and < 0.5 
M 5 . 

The transitional disk frequency for Coronet Cluster 
members may be spectral type dependent. At most 
one and probably none of the disks around stars ear- 
lier than K5 appear to be transitional disks (0-1/5). 
Stars between K5 and M2 also mostly have primordial 
disks (7/10) with 2-3 stars harboring transitional disks 
and 0-1 with debris disks. Accounting for the uncer- 
tain status of debris disk candidates, the relative fraction 
of transitional disks to all protoplanetary disks for stars 
M2 or earlier with plausible masses equal to 0.5 M Q or 
greater is then f(TD)/(f(TD)+f(PD)) » 0.15-0.27 (2/13, 
4/15). For K5-M2 stars only, 22-30% of the protoplan- 
etary disks are transitional disks. Transitional disks are 
more prevalent around stars later than M2, comprising 
50-58% of the protoplanetary disk population (5/10 or 
7/12). This high frequency is consiste nt with previous 
results from ISicilia-Aguilar et all (|2008f ). 

However, given the small population of Coronet Clus- 
ter stars studied here it is unclear whether our analy- 
sis truly identifies a statistically significant stellar-mass 
dependent frequency of transitional disks characteristic 
of most young clusters. Furthermore, it is possible that 
some transitional disks provisionally identified from mod- 
eling optical-to-mid IR data may instead be primordial 
disks whose weak emission is due to the growth and set - 
tling of submicron-sized grains (e.g. IFurlan et all [2006). 
Analyzing data from clusters where far-IR and submm 
data is available is required to more conclusively break 
these degeneracies. To explore these issues further, and 
provide a context for our results, we now analyze disks 
in other clusters. 

4. TRANSITIONAL DISK FREQUENCIES AND LIFETIMES 
AS A FUNCTION OF STELLAR MASS: A COMPARISON 
WITH TAURUS, IC 348, NGC 2362, AND n CHA 

4.1. Sample and Wavelength Ranges Used For Analysis 

5 In making these divisions, we are grouping together 
intermediate-mass (M* ~ 2-4 Mq) and solar-mass T Tauri stars. 
Grouping together these stars as the first spectral type bin may 
hide additional evolutionary trends since disk evolution for inter- 
mediate mass stars (> 2—3 Mq) and sol ar-mass stars may procee d 
at different rates even by 1-3 Myr (e.g. Currie and Kenyon 2009). 
Any differences between 1-4 Mq stars and lower-mass stars then 
applies only to the ensemble average of 1-4 Mq stars. 



To investigate the evolution of the transitional disk 
frequency with time, we analyze data in other clusters 
spanning an age range of 1-8 Myr. Comparing our results 
for the Coronet Cluster with 1-2 Myr-old Taurus and 
2-3 Myr-old IC 348 addresses whether a stellar mass- 
dependent transitional disk frequency is a general feature 
of 1-3 Myr-old clusters. By comparing these results with 
those for 5 Myr-old NGC 2362 and 6-8 Myr-old r) Cha, 
we investigate how the frequency of transitional disks 
evolves with time and estimate the typical transitional 
disk lifetime. 

Taurus - Members of the Taur us- Auriga star- f ormin g 
region are lis ted in catalo gs from IKenvon et all (|2008h . 
iRebull et all (|2010f ). and ILuhman et al.l (120091) which 
are updates of the classic IKenvon and Hartmannl (|1995[ ) 
catalog. We reanal yze Taurus members studied in 
ILuhman et al.l (|2010t ) to provide direct comparisons 
with their disk analysis. For Spitzer phot ometry, we use 
the IR AC and MIPS 24 \xm data from ILuhman et al.l 
(120101) and MIPS 70 fj,m and 160 \im data from 
Rebull et al.l (I2010D. W e adop t optical photomet r y from 



Kenyon and Hartmannl til995l) . IWhite and Ghezl (|2001l) . 
and [Audard c t al.l (120071 ) and submillimeter measure- 



ments f rom I Andrews and [Williams! (I2005D. Uewitt et~al 



(|1994[) . I Jensen et al.l (|1994D . and iBeckwith and Sargentl 
(1991). To dercddcn stars, we determine optical 



extinct ions for Ta urus memb ers, using est imate s 

from IFurlan et al.l (120061) 



lRebull et al.1 (I20T0I) 
IKenvon and Hartmannl ( 19951 ). and ILuhman et al.1 
(|2010T l as starting points and then fit the optical/near- 
IR portio n of the SEP to a synthetic SED produced 
using the ICurrie et al.l (1201 Ol) intrinsic colors and the 
iCurrie et al.l (|2010l ) and ILuhman et~aT1 (|2003l) effective 
temperature scales. 

A large number of Taurus sources have complete SEDs 
from optical to far-IR/submillimeter wavelengths. We 
model the SEDs of 25 K5-M6 stars whose IR colors plau- 
sibly identify both primordial disks and transitional disks 
based on var ious color-color selection criteria in the lit- 
eraturc (e.g. Lada ct al. 2006; Currie and Kenyon 2009; 
ILuhman et all 120101; rMuzerolle et al.ll2010l ) 6 . Quantita- 
tively, our selection criteria includes K5-M6 sources with 
K s -[8] (dcrcddcncd) = 1.25-2. We further restrict our 
sample to sources whose far-IR/submillimcter data can 
yield good constraints on the range of plausible disk 
masses (see Appendix). Specifically, we require that the 
source has published submillimeter data or a detection in 
the far-IR MIPS bandpasses (70 /im or 160 fim). While 
this criterion may bias our selection against including the 
most depleted disks, especially for the lowest mass stars, 
we include it to be conservative. 

IC 348 - For IC 348, we compile optical, near-IR 
and Spitzer/IRAC ph otometry, s p ectra l types and ex- 
tinctions presented in lLada et all (j2006L and references 
therein) a nd MIPS photomet r y and upper limits pre- 
sented in ICurrie and Kenvonl (|2009l ) for 307 members 
listed in both studies as well as data for 41 additional 
spectr oscopically confirmed members from lMuench et al.l 
(2007). We restric t our analysis to so u rces w ith MIPS 
24 \xm detections. ICurrie and Kenvonl ((2009) identified 
mid-IR colors for disks in different evolutionary states 

6 We do not present modeling results for sources such as DM Tau, 
whose classification as a transitional disk has broad agreement 



7 



by modeling the SEDs of newly-detected disks in IC 348, 
but their sample is comprised mainly of stars later than 
M3-M4. Here, we model the SEDs of 25 I C 348 stars pre- 
viousl y detected bv lLada et al.l (|2006T ) and lMuench et all 
(2007) with full SEDs through 24 /im focusing on those 
with the same range of mid-IR colors analyzed in Taurus. 
Most of these 25 stars are earlier than M3-M4. 

As with the Taurus sample, we use the published 
extinction estimates as a starting point and fit the 
optical/near- IR portion of the SED to arrive at a final 
value. In most cases, t he best-fit optical extinction (Ay) 
matches that listed bv lLada et al.1 (|2006l ) . Exceptions in- 
clude IDs 26, 58, 97, 110, and 9024: their new extinction 
estimates are generally larger. 

NGC 2362 - For NGC 2362, we analyz e the Spitzer 
IRAC a nd MIPS photomet ry presented by ICurrie et al.l 
(|2009t) . ICurrie et all £2009) present data for two mem- 
bership lists for NGC 2362 : members/probable members 
identified byjDahm (2005) and candidate members iden- 
tified by llrwin et al.l ((2008). Here we analyze Spitzer 
data only for members/probable members identified by 
iDahml (|2005l ) to be conservative. We include the two 
sources - IDs 41 and 63 - not considered in analysis by 
iLuhman et al.l (|2010f ) . because both have » 5c detections 
a nd thus should be in cluded in our analysis (Figure [9]). 

ICurrie et al.l (|2009h already model the SED s for NGC 
2362 stars with the iRobitaille et al.l ([2006) radiative 
transfer grid to identify disks with inner holes. Therefore, 
we only compare NGC 2362 SEDs through 24 fim to the 
appropriate flat reprocessing disk models. As starting 
points for our fla t disk compar i sons, we use extinction es- 
timates listed in ICurrie et al.l (l2009Tl that were based of f 
of comparisons with the iKenvon and Hartmannl (119951) 
colors and off of output from fits with the lRobitaille et all 
models treating the extinction as a free parameter. 

■q Cha - For 77 Cha, we use the IRAC and MIPS 
24 fi m and 70 iim pho t ometry and spectral types 
from ISicilia-Aguilar et al.l (|2009j i. MIPS 70 fim up- 
per l imits and MIPS-160 fim photometry/upper limits 
from [ Gautier et al. (200§), and IRS data presented i n 
iBouwman et al.l (|2006l ) and ISicilia-Aguilar et al.l (|2009l ). 
To model ij Cha sources, we compare source SEDs to 
appropriate flat disk models and fit their SEDs through 
160 fim with the Robitaille model grid to identify inner 
holes and estimate disk masses. 

4.2. Classification Method 

Our disk classification method is an expanded, more 
thorough version of that used in §3. As in §3, we iden- 
tify transitional disks with inner holes from fitting source 
SEDs with the Robitaille radiative transfer models. The 
difference between our method here and in §3 is in our dif- 
ferentiation between primordial disks and homologously 
depleted transitional disks. In §3, members of the lat- 
ter group were identified as those having optically thin 
mid-IR emission (from comparing their emission to the 
flat disk limit) but lacking evidence for inner holes. As 
stated in §3, we considered the classification of these 
disks as provisional since weak mid-IR emission could be 
due to a low overall mass of submicron to submillimctcr- 
sized dust or the growth and settling of just the smallest, 
submicron-sized grains to micron sizes or larger. Model- 
ing based on 1-24 fim data alone is unable to break this 
degeneracy. 



However, SED modeling of sources with far-IR/submm 
data helps break this degeneracy (e.g. iWood et al.| [2002') . 
Taurus and 77 Cha include many sources with far-IR 
and/or submm data. Thus, SED modeling of Taurus and 
77 Cha sources can identify which ones with optically- 
thin mid-IR emission but no inner hole have 1) a low 
inferred disk mass or 2) an inferred disk mass compa- 
rable to optically-thick primordial disks. The former 
group are classified as ho mologously depleted transitional 
disks ()Currie et al.ll2009h : the latter group could be pri- 
mordial disks un dergoing extreme gr ain growth/dust set- 
tling/shadowing (|Furlan et al.ll2006h . 

To identify the range of disk masses for optically-thick 
primordial disks and thus most appropriate limit for di- 
viding between primordial disks and homologously de- 
pleted transition al disks, we refer to SEP m odeling re- 
sults for Taurus ((Andrews and W illiams 2005) which are 
discussed in detail in the Appendix. Briefly, typical pri- 
mordial disk masses as a fraction of stellar mass are 
Mdisfc ~ 0.01 M* for Taurus: the interquartile range 
covers M disk ~ 0.003-0.03 M* (Appendix, Figure 7 
I Andrews and Williamsll2007ft 7 . About 92-96 % of the pri- 
mordial disks in the I Andrews and Williams! Taurus sam- 
ple have inferred masses greater than 0.001 M* (cf. Ap- 
pendix). As discussed in the Appendix, far-IR data (A 
> 70 fim) is sufficient to estimate the disk mass at a 
precision useful for our study. 

Thus, if we define the range of primordial disk masses 
by the interquartile range, the appropriate limit differ- 
entiating optically-thick primordial disks from homolo- 
gously depleted transitional disks is Mdisk — 0.003 M*. 
If we use the full range of primordial disk masses, the 
appropriate limit is Mdisk — 0.001 M*. As discussed 
in the Appendix, these estimates for typical masses 
of optically-thick primordial disks, and thus our crite- 
ria for separating primordial and transitional disks, are 
not compromised by incompleteness. In our analysis, 
the dis k mass is the median value of masses from the 
best-fit IRobitaille et al.l (|2006l) mo dels; to determ i ne the 
fractional disk mass, we use the iBaraffe et all (|1998f ) 
isochrones and th e effective temperatu re scales from 
ICurrie et all (|2010D for K5-M2 stars and ILuhman et~aT1 
((2003J) for later stars. 

Combining our criteria from §3 with additional criteria 
presented in this section, we consider a disk in Taurus or 
7/ Cha to be a primordial disk if it meets the following 
two conditions: 

• It lacks evidence for an inner hole/gap as inferred 
from SED modeling. 

• Its mid-IR emission exceeds the optically-thick, flat 
disk limit or it has an inferred disk mass of M<j > 
0.003 M*. 

Disks with inferred masses of greater than 0.001 M* and 
up to 0.003 M* but no inner hole are considered to be bor- 
derline cases, as their classification depends on whether 
the interquartile range or full range of optically-thick pri- 
mordial disk masses is used to separate primordial disks 

7 [Andrews and Williams ( 20051. l2007h use the term " Class 11 ob- 
jects" to describe sources with optically-thick disk emission, though 
their sample includes some transitional disks with inner holes (e.g. 
UX Tau, CoKu Tau/4, GM Aur) 



from homologously depleted transitional disks. Disks 
with masses < 0.001 M*, optically-thin mid-IR emission, 
and no inner hole are classified as homologously depleted 
transitional disks. We then determine the transitional 
disk frequency by dividing the number of transitional 
disks by the number of transitional + primordial disks. 
To analyze disks in clusters lacking far-IR/submm data 

- IC 348 and NGC 2362 - we first assign provisional 
classifications based on §3 methods, which yield provi- 
sional frequencies for transitional disks. Since mid-IR 
data alone is unable to differentiate between disks with 
a low dust mass and those with substantial dust settling, 
some the disks identified as transitional based on mid- 
IR data could instead be primordial disks. Thus, we use 
the Taurus results to identify the mid-IR colors of tran- 
sitional disks and determine the " contamination rate" of 
the transitional disk population selected by mid-IR data 
alone. Transitional disk frequencies for IC 348 and NGC 
2362 are then revised after taking into account this con- 
tamination rate. 

4.3. Transitional Disks in Taurus 

Modeling results for 25 Taurus members are summa- 
rized in Tables [5] and [7] Several of these members have 
disks with inner holes - e.g. UX Tau, LkCa 15, and 
J04333905+2227207. Many others - e.g. VY Tau - 
lack clear evidence for an inner hole but have emission 
falling below the flat, optically thick disk limit (Figure 
[T0|) . They also have weak submm emission (Table [7]). 
The weak near-IR through submm emission from these 
disks identifies a substantial depletion of the entire dust 
grain population from hot inner disk regions to cold outer 
disk regions. Nearly all sources (94 %) with K s -[8] < 1.75 
have emission weaker than the optically-thick, fiat disk 
limit. 

4.3.1. Best- Estimated Masses of Disks in Taurus 

Figure [TTJ displays the SEDs and histogram distribu- 
tion of masses from best-fit models for several disks that 
lack inner holes but have weak near-to-mid IR emission 

- VY Tau, JH-223, ZZ Tau, and FP Tau - and indi- 
cates that they have low inferred masses (Mdisk < 0.001 
M*). Such sources fit our criteria for being a homol- 
ogously depleted transitional disk. Several others (e.g. 
V807 Tau) lacking inner holes have optically-thin near- 
to-mid IR emission and Mdisk = 0.001-0.003 M*: these 
are borderline cases. The intrinsic sampling of the Ro- 
bitaille grid is heavily weighted towards disk models with 
masses ~ 0.001-0.01 M© (shaded region). In spite of this 
bias, the inferred disk masses for sources like those in 
Figure [TT] are at least an order of magnitude lower. 

Mass estimates based on submillimeter data alone 
(jAndrews and Williams 2005[) provide a simple but in- 
dependent check on the estimates derived from radia- 
tive transfer modeling. All of the homologously depleted 
transitional disks with submillimeter data that are listed 
in Table [7] have non detections at 350-850 /xm, and thus 
their masses have a wide range of possible values. How- 
ever, the submm-derived upper limits constrain the mass 
of optically thin emitting dust, which is useful since 
all of these sources have optically thin near-to-mid IR 
dust emission and optical depth decreases with increas- 
ing wavelength. For the homologously depleted disks, 
submm-derived upper limits range from 5.5 xl0~ 4 M* 



to 10~ 3 M*. In other words, submillimeter data alone 
provides evidence that these disks are less massive than 
primordial disks. 

4.3.2. Frequency of Transitional Disks in Taurus 

To compute the relative frequency of transitional disks 
around K5-M6 stars, we combine the trans itional disk 
population identified bv lLuhman et ail (|2010| ) and newly 
identified transitional disks from Table [7] For the pri- 
mordial disk population, we consider all Taurus members 
with Spitzer photometry that are not C l ass I protostars. 
Based on mid-IR colors, iLuhman et al.1 ([20101 ) identified 
eighteen transitional disks around K5-M6 stars: nine 
around K5-M2 stars and nine around later stars. To 
these, we add 6-7 transitional disks around K5-M2 stars 
and 3-7 transitional disks around later stars identified 
from SED modeling. 

Combining both transitional disk samples, the fre- 
quency of transitional disks is 0.19-0.22 (27-32/145) for 
the entire population, 0.20-0.21 (15-17/77) for K5-M2 
stars, and 0.18-0.24 (12-16/6 8) for M3-M6 stars Two 
transitional disks identified bv lLuhman et al.l ()2010l ) have 
extremely tenuous 24 fim excess emis sion (K,,-[24l ~ 
0.5-1 ) and could be young debris disks ([Carpenter et al.l 
2009). If we remove these sources, the transitional 
disk frequency for the entire population drops to 0.18 
and 0.22, while the frequency for K5-M2 stars drops 
to ~ 0.16-0.18. Thus, we find a slightly higher per- 
centage of transition al disks around K5-M6 stars than 
ILuhman et ail (|2010h find for K5-M5 stars (~ 20% vs. 
their 13% (15/113)). 

4.3.3. Mid-IR Colors of Primordial Disks and 
Transitional Disks 

To guide our analysis of other clusters lacking far- 
IR/submm data (e.g. IC 348, NGC 2362), we determine 
the mid-IR colors (K s -[8]) that identify primordial disks 
and transitional disks. Figure [T^] plots the mid-IR col- 
ors of Taurus members separated by disk evolutionary 
state. Only members with primordial disks (black cir- 
cles) have K s -[8] > 1.75. Taurus members with K s -[8] 
bluer than the horizontal dashed line (K s -[8] = 1.25) do 
not have primordial disks. Members with intermediate 
colors include primordial disks, transitional disks (grey 
circles enclosed by squares), and borderline cases (black 
circles enclosed by squares). Thus, there are two regions 
of mid-IR color-color space containing transitional disks: 
an " uncontaminated" region with blue colors consisting 
only of transitional disks and a "contaminated" region 
with redder colors which includes both transitional disks 
and primordial disks. 

Comparing the number of transitional disks and pri- 
mordial disks in the " contaminated" color region for Tau- 
rus members with far-IR/submm data provides an es- 
timate for the level of contamination in other clusters 
lacking far-IR/submm data. As shown by Figure [T2l and 
TableEl there are 9 K5-M2 members with K s -[8] = 1.25- 
1.75: 6 sources identified as transitional disks, two iden- 
tified as primordial disks, and one borderline case (V807 
Tau). Thus, 22-33% (2-3/9) of the sources with K s - 
[8]=1.25-1.75 and provisionally identified as transitional 
disks from mid-IR data alone could in fact be primordial 
disks with substantial dust settling. For M3-M6 stars 
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with K a -[8]=1.25 -L75 3 3 sources liste d in Table plus 
two identified bv lLuhman et al.1 (|2010f ) - J04330945 and 
J04213459 - are transitional disks, one is a primordial 
disk, and four are borderline cases. Thus, the contami- 
nation rate for M3-M6 stars with K s -[8] = 1.25-1.75 is 
10-50% (1-5/10). 

4.3.4. Incorporating Uncertainties in Disk Mass 
Estimates 

Here we quantitatively assess how the intrinsic un- 
certainty in disk mass estimates affects our conclusions 
about the frequency of transitional disks in Taurus and 
the contamination level from selecting transitional disks 
by mid-IR colors. First, we compare the interquartile 
range of disk masses from best-fit models for each star. 
Figure [13] shows that the range of disk masses for homol- 
ogously depleted disks and primordial disks are distinct 
(see also Table [S]). The upper quartile of disk masses for 
homologously depleted disks ranges from ~ 9xl0 -4 M* 
to 4xl0 -3 M*. The lower quartilcs of masses for pri- 
mordial disks is systematically larger and does not over- 
lap, ranging from 4xl0 -3 to 10~ 2 M*. The masses 
of disks in borderline cases are also systematically low, 
though the upper quartiles for three of the five stars - 
J04231822+264116, J04153916+2818586, and ITG 15 - 
overlap slightly with the lower quartile masses for several 
primordial disks (e.g. GH Tau). 

Second, we consider the full distribution of disk masses 
from best-fit models for each star to identify the frac- 
tion of models yielding M disfc /M* < 0.001 and 0.003 (f 
< 1, 3xl0~ 3 ), our thresholds identifying homologously 
depleted transitional disks and borderline cases. Our 
goal here is to derive a purely probabilistic estimate of 
the transitional disk population from T able [6] to com- 
bine w ith the transitional disks listed by iLuhman et al.l 
(2010) and then to compare with previous estimates of 
the total population in this section. Here, the number of 
transitional disks around n stars listed in Table [6] each 
with f(TDj of their disk models fulfilling our definition 
of a transitional disk is 

n 

n(rr>) = 5^/(TA). (i) 

For disks with inner holes - e.g. LkCa 15 - determined 
by imaging and SED modeling f(TD^) =1, otherwise it 
equals f(< 1, 3 x 10~ 3 ) as listed in Table[Sl This exercise 
yields n(TD) = 10.6-13.8 or ~ 11-14 (n(TD) = 5.5-7.1 
for K5-M2 stars, 5.2-6.7 for M3-M6 stars). 

Combining our probabalistic estimate for the tran- 
sitional disks population drawn f rom Table EH and the 
18 transitional disks identified by ILuhman et al.l (|2010f ) 
yields total population of 29-32 (18+ [11,14]), in agree- 
ment with estimates derived in §4.3.2. The estimates 
divided by spectral type range are likewise similar: 15- 
16 (9+[6,7]) for K5-M2 stars and 14-16 (9+[5,7]) for 
M3-M6 stars. Finally, the resulting contamination esti- 
mates agree well with previous ones: 22-33% (2-3/9) for 
K5-M2 stars and 10-30% (1-3/10) for M3-M6 stars 8 . 

8 For K5-M2 stars, there are nine sources within our contamina- 
tion region, our exercise predicts that ~ 6—7 are transitional disks, 
so the contamination rate is 2—3/9. For M3— M6 stars, there are 10 
sources within the contamination region. Two of these are identi- 



4.4. Frequency of Transitional Disks in IC 34-8 

Table [5] summarizes our modeling results for selected 
IC 348 members; SEDs for members representing a range 
of mid-IR disk emission and provisional evolutionary 
states are shown in Figure [T4l As with Taurus, we find 
that all but one of the members with K s -[8] < 1.75 have 
optically-thin mid-IR disk emission. Since our Taurus 
analysis finds that most members with these colors have 
low disk masses or inner holes (29-32/35), we infer that 
most IC 348 members with these colors also have low disk 
masses or inner holes. Of the 25 IC 348 stars modeled, 
17 are consistent with being transitional disks based on 
1-24 fim data, while 8 are primordial disks. 

To determine the frequency of transitional disks in IC 
348, we use contamination estimates based on mid-IR 
colors from Taurus. Following the Taurus results, we 
adopt K s -[8] = 1.75 as our fiducial division between pri- 
mordial disks and transitional disks and derive a provi- 
sional frequency for transitional disks. Then we derive a 
revised frequency taking into account contamination by 
primordial disks with K s -[8] = 1.25-1.75. 

Based on this method, 8/33 K5-M2 stars and 43/81 
M2.5-M6 stars have evidence for a disk with K s -[8] < 
1.75, and thus fit our nominal near-to-mid IR transitional 
disk criteria. Of the 8 (44) disks around K5-M2.5 (M2.5- 
M6) stars provisionally classified as transitional, 4 (19) 
have mid-IR colors within the "contaminated" region. 
Thus, corrected for contamination, about 32-40% of K5- 
M6 IC 348 disks are transitional disks. The transitional 
disk population is heavily weighted towards later-type 
stars: 18-24% of disks around K5-M2 stars have transi- 
tional disks (6/33 to 8/33), while ~ 42-51% of M2.5-M6 
stars have transitional disks (34/81 to 41/81). 

4.5. Frequency of Transitional Disks in NGC 2362 

Table [10] summarizes our modeling results for NGC 
2362. Near-to-mid IR emission from most cluster mem- 
bers lies well below the optically thick, flat reprocessing 
disk limit, indicating that the mid-IR emission is at least 
marginally optically thin (r < 1). Comparing this mod- 
eling to radiative transfer m odeling performed for NGC 
2362 in ICurrie et al.l (|2009h and for Taurus in this work 
shows that many sources with weak/negligible 1-8 \xm 
emission and optically-thick 24 fim emission require inner 
holes while many with marginally optically-thin near-to- 
mid IR emission have low inferred disk masses. Figure [TBI 
displays 4 representative SEDs comparing the observed 
emission to the optically thick, flat disk limit; 

Our results are consistent with those from lCurrie et all 
(2009) who argue that transitional disks ar e equal to or 
greate r in number than primordial disks. ICurrie et al.l 
(2009) find that 81% of disks around K0-M3 members 
are transitional disks (17/21). Here four disks are pri- 
mordial disks (IDs 111, 139, 187, and 202), and one (ID 
177) appears to be a borderline case (since it may have 
an inner hole) but is counted as a primordial disk in 
the analysis below to be conservative (since it has K s - 
[8] > 1.75). The rest (16; 14 of which are K5-M3 stars) 
have near-to-mid IR emission consistent with being tran- 
sitional disks. Of these, five K5-M3 stars either have 

fied as transitional disks from Luhman et al. (2010), whereas our 
exercise predicts that 5—7 out of the eight others are transitional 
disks. Thus, the contamination rate is 1—3/10 
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K s -[8] < 1.25 (thus lying outside the range of colors con- 
taminated by primordial disks) or show evidence for in- 
ner holes. The other 9 K5-M3 stars have mid-IR col- 
ors between K s -[8] = 1.25 and 1.75, within the contam- 
inated region. Based on our Taurus analysis, 2 K5-M2 
stars and 0-1 M2.5-M3 stars provisionally identified as 
having transitional disks instead have primordial disks 
with dust settling. Therefore, we arrive at a final tran- 
sitional disk frequency of 58%-63% (11/19-12/19) for 
K5-M3 stars. The frequency for K5-M2 stars (63%) and 
M2.5-M3 stars (33-66%) are statistically indistinguish- 
able given the small number in the latter group (3). 

4.6. Frequency of Transitional Disks in r\ Cha 

Table [Til summmarizes our modeling results for rj Cha. 
Of the 15 late-type members with IRAC and MIPS pho- 
tometry, 8 have evidence for a disk. Only three sources 
(RECX-11, J0843, and J0844) fulfill our definition of 
a primordial disk. Three members have transitional 
disks with inner holes (RECX-3, 4, and 5); two have 
homologously depleted transitional disks (RECX-9 and 
J0841). Thus, rj Cha has a high transitional disk fre- 
qu ency (5/8 or 63%). Our results agree with those 
of ISicilia-Aguilar et al.l (|2009| ) who separately analyze 
Spitzer IRS data to find that transitional disks comprise 
50-75% of the disk population. A "probablistic" esti- 
mate of the transitional disk fraction (e.g. as in Section 
4.3) yields the same answer: n(TD) — 4.95 or ~ 5. 

Gas diagnostics are consistent with our argument that 
most disks in rj Cha are transitional disks. RECX- 
5 and 9 are accreti ng at extremely low rates (4-5 
x 10- 11 Mq yr" 1 ) (TLawson et all l200l . RECX-3, 
RECX-4, and J0841 lack evidence for accretion cir- 
cumstellar gas. RECX-3, 5, and 9 lack evidence for 
rovibrational H2 emission from wa rm circumstellar gas 
(jRamsav-Howat and Greaves! 120071 ) . Only J0843, iden- 
tified as a primordial disk from SED modeling, shows 
evidence for H2 emission and for substantial accretion 
(10- 9 M yr" 1 ). 

4.7. Revised Transitional Disk Frequency for the 
Coronet Cluster 

Based on our analysis of Taurus data in §4, some frac- 
tion of Coronet Cluster stars determined to have transi- 
tional disks based on mid-IR data probably are primor- 
dial disks with dust settling. In particular, G-l (MO), 
G-14 (M4.5), and G-87 (M1.5) have K s -[8] colors placing 
them within the contaminated region: other transitional 
disks have inner holes or lie within the uncontaminated 
region. Accounting for contamination, the size of the 
transitional disk population drops to 11-22% for K5-M2 
stars (1/9, 2/9) and 30-50% for M2.5-M6 stars (3/10, 
6/12), where the ranges account both for contamination 
by primordial disks with substantial dust settling and 
for possible misclassification of debris disk candidates. 
Thus, our transitional disk frequencies for all five clus- 
ters now account for uncertainties in disk classification 
based on mid-IR data alone. 

4.8. The Evolution of the Transitional Disk Population 

Combining our results for the Coronet Cluster, Tau- 
rus, IC 348, NGC 2362, and rj Cha, we now investi- 
gate the frequency of transitional disks as a function of 



age and stellar properties. Table [12] lists the relative 
frequency of transitional disks as a function of time in 
three spectral type bins: mid-K to mid-M stars (K5- 
M6), K5-M2 stars, and M3-M6 stars. To compare 
with theoretical predictions, we overplot the approxi- 
mate loci (Figure IT7l top panel, solid and dotted lines) 
expected from the lAlexander and Armitagel (|2009D disk 
evolution/ planet formation models, which yield frequen- 
cies for transitional disks as a function of time and as- 
sume a disk clearing timescale of ~ 0.5 Myr. 

The transitional disk frequency for K5-M6 stars in- 
creases as a function of time from ~ 20% at 1- 
2 Myr to ~ 50-60% by 5 Myr. The frequencies 
are consist ently 1.5-3 times highe r than those pre- 
dicted by lAlexander and Armitagel (|2009f ) for a 0.5 
Myr transition timescale . Wi thin the context of the 
lAlexander and Armitagel (|2009f ) study, our results indi- 
cate that the transitional disk phase on average lasts 
longer than 0.5 Myr. 

The bottom panel of Figure 1171 compares o ur frequen- 
cies with those from IMuzerolle et al.l (|2010l ) who use a 
more empirically-based approach to identify transitional 
disks based primarily on IRAC and MIPS flux slopes 
(bottom panel). Our transition al disks frequenc ies are 
substantially larger than those IMuzerolle et al.l derive 
based on their self-described classical definition for what 
constitutes a transitional disk (diamonds): a disk with 
an inner hole a nd optically-thick outer disk. However, 
IMuzerolle et alJ argues for an expanded definition for a 
transitional disk, including disks whose mid-IR emission 
plausibly identifies a more homologous depletion of emit- 
ting dust and disks with inner holes/cavities and more 
optically thin outer disks. Adopting this definition, their 
transitional disk frequencies (triangles) are in good agree- 
ment with those we derive for clusters older than 2 Myr. 
For the two clusters we both analyze, IC 348 and rj Cha, 
our frequencies are nearly identical. 

To derive a timescale for the transitional disk phase, 
we follow a parametric, Monte Ca rlo approach, similar to 
that from Muzcroll e et al.l l)2010l ). evolving a population 
of 10 5 stars to simulate the frequency of disks in differ- 
ent states for clusters with ages less than 10 Myr. We 
adopt a 3 Myr e-folding timescale for the protoplanetary 
disk stage, consistent with Spitzer observations (Figure 
[T51 top panel). Some fraction of the protoplanetary disk 
lifetime is spent in the transitional disk phase. To assess 
the duration of this phase, we vary the transitional disk 
e-folding timescale (trc); correspondingly, the e-folding 
timescale for the primordial disk phase is 3 Myr - Itd- 
Through the first 0.54 Myr, we fix the disk fracti on at 1 
to ac count for the lifetime of Class I protostars (jEvansI 
12009t h Finally, we assume a 1 Myr age dispersion at each 
cluster age. 

As shown by Figure [18] (bottom panel) , transitional 
disk frequencies in all clusters are most consistent with a 

1 Myr transition timescale, lying well above predictions 
for a 0.1-0.5 Myr timescale and below predictions for a 

2 Myr timescale. Thus, our analysis confir ms earlier re- 
sults b ased on NGC 2362 and rj Cha fr om ICurrie et al.l 
l|2009h and ISicilia-Aguilar et all (|2009l ) that the mean 
transitional disk lifetime must be an appreciable frac- 
tion of the total p rotoplanetary disk lifetime, estimated 
to be ~ 3-5 My r dHernandez et all l2007at ICurrie et all 
2009). Adopting IMuzerolle et all (|2010f )'s model for the 



11 



frequency of transitional disks with time yields the same 
qualitative result. They predict smaller transitional disk 
frequencies at 3-10 Myr than we do for a 1 Myr tran- 
sition timescale. The identification of a long transition 
timescale based on our predictions is then conservative. 
On the other hand, our results are in conflict with other 
recent claims, particularly those of lLuhman et al.l (|2010f ) 
who argue for a transitional disk lifetime less than 0.5 
Myr. In §5, we compare our methodology with their's 
and those from other recent studies. 

We emphasize that our disk frequencies account for 
disks with both transitional disk morphologies: disks with 
inner holes and homologously depleted disks. Homolo- 
gously depleted transitional disks generally appear to be 
more numerous than those with inner holes, up to a fac- 
tor of 1.5-2 larger in number for the Coronet Cluster and 
NGC 2362. However, these differing frequencies do not 
constrain the relative lifetimes of disks following these 
evolutionary paths: we cannot tell between a morphol- 
ogy with a shorter lifetime and one that operates less 
frequently. We also note that these results are statisti- 
cal. In other words, while there is some uncertainty in 
the state for individual objects, these uncertainties are 
averaged out at the end. 

At least for the Coronet Cluster and IC 348, the transi- 
tional disk frequency is spectral type dependent (Figure 
IT5|) . The transitional disk frequencies for M3-M6 stars in 
these clusters is marginally but consistently higher than 
those for K5-M2 stars by about a factor of ~ 1.5. For- 
mally, M3-M6 stars in Taurus also more frequently have 
transitional disks than do K5-M2 stars, though the dif- 
ferenc e is not statistically significant. iMuzerolle et al.l 
(2010) also find that the transitional disk population is 
dominated by later-type stars, though their correlation 
is stronger than what we find. 

A spectral-type dependent transitional disk frequency 
may help e xplain the high tran si tional disk frequency 
found by ISicilia-Aguilar et al.l (|2008l) . Since the 
ISicilia-Aguilar et al.l ( 20081 ) sample is drawn heavily from 
stars later than M2, their sample is predisposed to- 
wards having a high transitional disk fraction compared 
to cluster-averaged values and compared to those for 
brighter, more easily detectable K to M2 stars. The tran- 
sitional disk frequency for the bulk cluster population 
(e.g. including members in Table [3]) is similar to that for 
other 1-3 Myr-old clusters. At least one other cluster, 
IC 348, also has a very high transitional disk fraction for 
very low-mass stars. 

5. DISCUSSION 

5.1. Summary of Results 

This paper expands upon previous Spitzer analysis of 
the Coronet Cluster disk population around v er y low - 
mass stars performed by ISicilia-Aguilar et al.l ((20081 ). 
presenting new photometry for other cluster members, 
including many with higher stellar masses. We use 
SED modeling and comparisons with simple disk mod- 
els/empirical metrics to assess the Coronet Cluster disk 
population, focusing on identifying and characterizing 
candidate transitional disks. By analyzing the disk pop- 
ulation for other 1-8 Myr-old clusters, we investigate 
the utility of using mid-IR data to probe disk evolution- 
ary states, determine how the frequencies of transitional 
disks change with time, and investigate how this lifetime 



may depend on stellar properties. Our study yields the 
following major results: 

• The Coronet Cluster contains a high frequency 
of transitional disks: ~ 30% (formally, 21-36%) 
of the disk population around K5-M6 stars. Al- 
though the samples are small, our analysis hints 
at a spectral type- dependent frequency for tran- 
sitional disks, confi r ming e arlier suggestions by 
ISicilia-Aguilar et al.1 (|2008l l2009| ): they appear 
more frequently arou nd the very low-mass star s 
that were the focus of ISicilia-Aguilar et al.l ((2008). 
IC 348 exhibits the same spectral type/stellar- mass 
dependent transitional disk frequency. 

• Based on optical to submm SED analysis of Taurus 
sources with high-quality photometry, we confirm 
that many disks lacking clear evidence for inner 
holes also have a low dust mass. Their mid-IR 
emission is due to a depletion of dust, not dust 
settling: these represent a more homologous deple- 
tion of disk material with time instead of an inside- 
out dispersal. Many transitional disks can only 
be identified from SED modeling as their mid-IR 
colors overlap with primordial disks. About 22- 
33% (10-50%) of disks around K5-M2 (M3-M6) 
stars identified as homologously depleted transi- 
tional disks based on optical to mid-IR data alone 
may in fact be primordial disks with a heavy deple- 
tion of submicron-sized dust but a large total dust 
mass. Disks with very different morphologies can 
occupy the same region of color-color space: SED 
modeling is required to avoid misclassifying disks. 

• Combining analysis of the Coronet Cluster disk 
population with that for Taurus, IC 348, NGC 
2362 and n Cha shows that the relative frequency 
of transitional disks around K5-M6 stars increases 
from - 20% at 1-2 Myr to - 50-60% by 5-8 Myr. 
Parametric modeling shows that this trend is im- 
plies a mean transitional disk lifetime near 1 Myr, 
not ~ 0.1-0.5 Myr. 

5.2. Transitional Disk Identification in Other Recent 
Work: The Importance of SED Modeling 

Other recent studies of transitional disks provide ev- 
idence that the weak-IR-to-submm dust emission from 
homologously depleted transitional disks is l ikely due to 
disk c learing, consistent with our results. ICieza et "all 
(2008) selected transitional disks as those with mid-IR 
emission weaker than the lower- quart ile Taurus SED and 
low levels of accretion, including many that lack any ev- 
idence for inner holes. They found that disks with weak 
mid-IR emission h ave low submmm-inferred disk masses. 
ICieza et al.1 (|2008[ ) interpret this trend as evidence in fa- 
vor of UV photoevaporation models for disk clearing: 
photoevaporation should begin to clear the inner disk 
once the total disk mass (probed by submm data) drops 
significantly. More generally, though, this trend implies 
that for some disks the mass of emitting dust drops si- 
multaneously over a wide range of stellocentric distance, 
c onsistent with a ho mologous depletion of disk material. 

ICieza et al.l (|2010f ) found that the disk mass and accre- 
tion rate for transitional disks, including those we would 
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identi fy as homologous ly depleted, appear to be corre- 
lated (|Cieza et al.ll2010l ) : disks with lower accretion rates 
or small upper limits have lower inferred disk masses. 
There is no clear reason why the accretion rate should 
drop because submicron-sized grains had grown; how- 
ever, disks undergoing a homologous depletion of their 
total mass in gas and dust should have low inferred 
masses and a lower frequency and rate of accretion. 

Other studies, partic ularly the C2D sur vey of transi- 
tional disks presented in lMerin et al.l (|2010t ). conclusively 
show that SED modeling is needed to properly charac- 
terize disks. By modeling optical through submm pho- 
tometry and spectroscopy of a mid-IR selected sample of 
candidate transitional disks ("cold disks" in their termi- 
nology), they determine which disks have inner holes and 
which lack evidence for an inner hole and investigate typ- 
ical mid-IR colors for both kinds of disks. Their analysis 
shows that the mid-IR colors of transitional disks overlap 
with primordial disks (e.g. their Figure 15). Computing 
the frequency of transitional disks using mid-IR colors is 
then error prone unless SED modeling of well character- 
ized sources is used to identify uncontaminated regions 
of color-color space and quantiatively assess contamina- 
tion where transitional disks and primordial disks have 
the same colors. This finding is qualitatively consistent 
w ith our Taurus model ing results. 

iLuhrnan et al.l (|2010l ) compile data for numerous 1-10 
Myr-old clusters to uniformly identify transitional disks 
and quantify their population self consistently. Their 
modeling follows a different approach than that followed 
here and in some previous work (e.g. ICurrie et alj [2009; 
iMerin et al.l l2010h . as they selec t synthetic K s - [ 8] and 
K s - [24] colors produced from the ID'Alessio et al.l (|2006f ) 
1+1D disk code as fiducial colors separating primordial 
disks from transitional disks. The model they use has 
two grain populations: where the "big" grain popula- 
tion extends from submicron to millimeter sizes and is 
confined to the disk midplane, while the "small" grains 
lie above the midplane and are < 1 /im in size. It as- 
sumes a negligible accretion rate (M = 1CP 10 M Q yr _1 ) 
and a "depletion factor", e=0.001, which removes 99.9% 
of the small dust grains. They find lower transitional 
disk frequencies and thus derive a shorter lifetime for 
the transi tional disk phase. They claim that previous 
work (e.g. ICurrie et all 120091 : ISicilia-Aguilar et al.l 120091 ) 
arguing in favor of a longer transitional disk lifetime im- 
properly identifies sources with emission weaker than the 
median Taurus SED as transitional disks or conflates a 
lack of near-IR e mission with evidenc e for disk clearing 9 . 

We agree with lLuhman et a l. (2010) that some primor- 
dial disks can have very blue colors due to dust settling. 
However, they did not quantitatively consider the pos- 
sibility that some transitional disks may (in some color- 
color diagrams) have redder colors than the bluest pri- 
mordial disks. Based on a more detailed, multiwave- 
length analysis, we derive higher transitional disk fre- 
quencies than they do for nearly all clusters, especially 
for those older than Taurus. Our disagreement primar - 
ily arises because many disks that ILuhrnan et al.1 ()2010P i 

9 A s can be seen from a simple inspection of the ICurrie et all 
(2009) analysis, their disk classification was instead based on com- 
paring SEDs t o the low er- quartil e Tau rus SED. This was clarified 
in later work ijCurrie and K cnvon 20091). 



identifies as being optically-thick primordial disks from 
their color criteria are consistent with having optically 
thin/marginally optically-thin emitting dust in the mid- 
IR and low inferred disk masses. If their weak mid-IR 
emission were due only to the growth and settling of the 
smallest grains, their far-IR to submm emission (and thus 
inferred disk masses) should be comparable to known pri- 
mordial disks, which is not observed. Therefore, these 
disks are more consistent with being transitio nal disks. 

There are other aspects of ILuhrnan et al.l 's analysis 
that weaken or under mine their conclusions. First, while 
ILuhrnan et al.l (|2010h 's adopted model results from a so- 
phisticated code, their conclusions are prone to large un- 
certainties, mostly because they use color comparisons 
from a single model to identify transitional disks, not 
full S ED modeling . Within the context of their model, 
the lLuhman et aTJ conclusions are extremely sensitive to 
assumed values for their free parameters, e and M, and 
thus highly selective. Assuming cither that disks remove 
99% of their small dust instead of 99.9% or have a small 
accretion rate of M = 10 -9 M Q yr _1 revises their fidu- 
cial mid-IR colors redwards by ~ 0.5 mags (e.g. K s -[8] 
- 1.75, K s -[24] - 4.5 for K5-M2 stars). Using these 
colors to classify di sks yields results sim ilar to those pre- 
sented here and in ICurrie et al.l (|2009l ) for 15 of the 16 
model combina tions of e (0.001 - 1) an d M (lO -10 -^ -7 
M Q ) studied in ID'Alessio et all (|2006l e.g. their Figure 
13). 

Second, while e may be a true free parameter, mea- 
sured accretion rates for many sources a re much higher 
than that from the disk model used by ILuhrnan et al.l 
(2010), rendering their adopted model inapplicable. For 
example, ID-8 5 in NGC 2 3 62, id entified as a transi- 
tional disk by ICurrie et alj (|2009h . is a likely accretor 
based on i ts H a emission. Using re l ations from [Dahm 
([2li0l and[H crczcg and Hillenbrand] ()2008f ) 10 . its accre- 
tion rate is ~ 4-7 x 10~ 9 M P) yr " 1 , or 40 to 70 times 
greater than what ILuhrnan et al.1 assu me. From com- 
paring the SED of ID-85 to output from ID'Alessio et al.l 
disk models with realistic accretion rates, this source 
sh ould be classifi ed as a transitional disk, contrary 
to ILuhrnan et all s results. Ta urus-Auriga, Tr umplcr 
37, and other clusters analyzed bv lLuhman et all include 
many sources with weak mid-IR emission (transitional 
disks / borderline cases ) whose measured accretion rates 
render ILuhrnan et al.l 's analy sis inapplicable, such a s ZZ 
Tau (1.3 x IP" 9 yr' 1 . iWhite and Ched [200ll see 
also ISicilia-Aguilar et al.ll2010D. 

Third, the model that ILuhrnan et "ail (|2010f l adopt 
likely 1) has an extremely low disk mass which 2) is al- 
most never optically thick and thus cannot distinguish 
optically-thick primordial disks from other dis ks. As 
sho wn in Tables 2.2 and 2.4 of lEspaillatl (|2009| ). where 
the ID'Alessio et al.l input parameters are detailed, the 
disk mass for models with parameters adopted by Luh- 
man is ~ 1.38-1.45 xl0~ 4 M : an order-of-magnitude 
(or more) smaller than our primordial/transitional disk 

10 Her e we assume E(B-V) = 0.01 (l/10th the cluster average 
listed by Moit mho et al.l r2001). an age of 5 Myr, a stellar radius 
for an M2 star of this age from Baraffc ct al. (1998), and Ri n = 
5 R *. The intrin sic uncertainty in the accretion rate is < 0.5 dex 
(sec Dahm 2008). According to SED modeling presented here, 
the reddening is probably higher, meaning that the intrinsic H Q 
luminosity is higher and the resulting accretion rate is also higher. 
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division for all but the latest-type stars and the me dian 
primordial disk mass ([Andrews and Williams! 120051 and 
Appendix). Furthermore, Figure 2.6 of |Espaillat (2009) 
shows that almost none of this disk is optically thick to 
its own r adiation (tr >> 1), unlike our adopted model 
and other iD'Alcssi oet al.l models. Thus, a major reason 
why the Luhman's adopted model yields weak mid-IR 
emission is because of its low dust mass, not just be- 
cause of d ust settling. To better parse disks into differ- 
ent states, iLuhman et al.l should have instead adopted a 
model yielding an optically-thick disk with a dust mass 
comparable to ty pical primordial disk masses (e.g. the 
fiducial model in lD'Alessio et aLll2006l except with e = 
0.001 and the accretion rate set to a negligible value). 

Finally, SED mode ling and imag ing of individual 
sources shows that the ILuhman et al.l color criteria fails 
to uniquely identify primordial/transitional disks. For 
example, UX Tau and Lk Ca 15 are iden tified as pri- 
mordial disks according to iLuhman et alJ 's color crite- 
ria as they are redder than the authors' primordial disk 
limit. But SED modeling and/or high-contrast imag- 
ing show that both have large gaps with optically-thin 
dust/no dust extending from 01 AU to ~ 50 AU 
from the star (iThalmann et al.ll20fol:lMulders et al.ll2010l; 
lEspaillat et al.l I2007D . SED modeling also shows that 
sources in NGC 2362 (e.g. ID-3), IC 348 (e.g. ID- 
97), and Trum pler 37 (e.g. ID 14-11) classified by 
ILuhman et al.l as primordial disks probabl y have inner 
holes and thus are transitional disks (e.g. ICurrie et al.1 
120091: ISicilia-Aguilar et al.ll2007l this work). 

Even though we disagree with ILuhman et alJ 's conclu- 
sions for reasons described above, their study reinforces 
a major point made here - mid-IR colors of disks can 
be highly degenerate. The degeneracy highlighted from 
both of our studies shows that full SED modeling is 
needed to accurately assess the frequency of transitional 
disks and the duration of the transitional disk phase. 
While time intensive, this approach provides de facto 
more information about disk properties than can be pro- 
vided by comparing source colors with synthetic colors 
resulting from a single disk model in a single color-color 
diagram. 

5.3. Transitional Disks In Young Clusters: Why the 
Paucity of Transitional Disks in Taurus Does Not 
Imply a Short Timescale for Disk Clearing 

Many pre-Spitzer studies have argued that the small 
number of transitional disks compared to optically-thick 
primordial disks and disk-less stars in young clusters 
like Taurus provides evidence in favor of a short (< 
0.1 Myr) transition disk phase due to rapid disk clear- 
ing (e.g. ISkrutskie et all Il990t ISimon and Pratol 119951 
IWolk and Walterlll996l and later work). Though transi- 
tional disks are infrequent compared to primordial disks 
and diskless stars in the youngest clusters (e.g. 20% for 
Taurus), it does not follow that the transitional phase 
must be very shortlived (e.g. < 0.1 Myr). Here, we re- 
view why this is the case, describe how our work and 
other recent studies support this position, and consider 
recent work tha t revives older claim s. 

As argued bv ICurrie" et~"aLl IpOOfll) . stars in Taurus are 
typically too young compared to the mean protoplane- 
tary disk lifetime (3-5 Myr) for the majority of them to 
begin clearing. Thus, it is unsurprising that relatively 



few show evidence of active clearing. Our parametric 
modeling quantitatively supports this argument, show- 
ing that transitional disks do not comprise the majority 
of disks in the youngest clusters even for lengthy transi- 
tional disk phases (e.g. Itd = 1 Myr). Thus, the size 
of the transitional disk population at ~ 1-2 Myr does 
not, by itself, clarify how long protoplanetary disks show 
signs of active clearing over a typical lifetime of 3-5 Myr. 
To explain the large population of d iskless stars in Tau- 
rus (~ 30% of the total population). ICurrie et all (|2009h 
argue that stars with close binary companions lose their 
protoplanetary disks far more rapidly than stars witho ut 
cl ose binary com panions (e.g. Ilreland and Krausll2008f ). 

ILuhman et al.1 disagree with both of these claims. 
They argue that age cannot explain Taurus 's paucity 
of transitional disks, since since diskless stars in Tau- 
rus are more dispersed and are thus olde r than stars 
with disks. We agree with ILuhman et al.1 (|2010T ) that 
an age spread may also explain the large number of disk- 
less sources in Taurus. However, this fact would simply 
confir m that that older stars are less likely to have disks 
(e.g. ISicilia-Aguilar et apl2006bHHernandez et aLll2007at 
ICurrie et al.ll2007l 120091 ) and thus does not provide evi- 
d ence in favor of rapid disk clearing. 

ILuhman et al.1 also argue that binarity cannot explain 
the large diskless population in Taurus. They note that 
the Taurus binary frequency declines with mass. They 
claim that if many Taurus disks are cleared by binaries 
then the disk fraction should be highest at low masses, 
which is not observed. However, the drop in binary fre- 
quency for low-mass stars is heavily weighted towards 
a drop in wide binary companions with sep arations of ~ 
200 AU or g reater, which are irrelevant here (|Kraus et al.1 
120061 120111 ). Moreover, many diskless stars in Taurus 
are tight binaries (Kraus et al. 2011, ApJ submitted). 
With tight binaries remove d, the popu l ation of diskless 
stars shrinks significantly. ICieza et al.1 (2009) also find 
that disks in short-period binary systems are less fre- 
quent than those around wide-separatio n binaries or sin- 
gle stars (see also iBouwman et "al"1l2006| ) . 

In summary, either binarity or an age spread explain 
the smaller number of Taurus stars with transitional 
disks relative to primordial disks and diskless stars. 
Moreover, parametric modeling shows that the frequency 
of protoplanetary disks in the transitional disk phase is 
perfectly consistent with an extended transitional disk 
phase. Thus, Taurus data provides no evidence that the 
transitional disk phase must be rapid. 

5.4. Limitations and Uncertainties of Our Study and 
Future Work 

5.4.1. Qualifications on Our Identification of 
Transitional Disks 

Our identification of some transitional disks (those la- 
beled as homologously depleted) is predicated on 1) our 
assumptions about the range of primordial disk masses, 
2) the reliability of the Robitaille models in estimating 
disk masses, and 3) the reliability of our assumed stel- 
lar masses. Current data shows that the lower quar- 
tile of masses for optically-thick primordial disks is 0.003 
M* and nearly all primordial disks have masses greater 
than 0.001 M*. Since not all Taurus members have been 
targed by submillimeter observations, it is possible that 
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future surveys will uncover more disks with optically- 
thick near-to-mid IR emission but inferred masses that 
are systematically lower than the interquartile range/full 
range from current observations. 

Additionally, as shown in Figure [HJ the best-fit disk 
models for homologously depleted transitional disks have 
median values for absolute disk masses that yield frac- 
tional disk masses lying below our adopted primordial 
disk limit. However, some of the models for each disk 
(typically, 10-30%) yield masses lying above 0.001 M*. 
Thus, it is possible, albeit unlikely, that some of these 
disks are not low mass. 

However, submillimcter data constrain the mass of 
emitting, optically thin dust to be less than these lim- 
its. Moreover, the best-fit models for some primordial 
disks and all borderline cases include those with masses 
below 0.001 M*. Thus, uncertainties in determining disk 
masses do not necessarily bias our results in favor of find- 
ing a high frequency of transitional disks. Similarly, the 
classification of some sources with optically-thin near-to- 
mid IR emission but large inferred disk masses (e.g. V836 
Tau) as primordial disks is conservative. While these 
sources lack clear evidence for cleaned inner holes/gaps, 
their inner regions may be optically thin, and they may 
be undergoing the early stages of inner disk clearing. 

More genera lly, the intrinsic sampling of the 
iRobitaille et al.l grid is non- uniform across all fitting pa- 
rameters, which in principle could undermine our con- 
clusions. As mentioned previously, disks without inner 
holes are more frequently represented than those with- 
out holes and massive disks (Mdisk ~ 0.01 M ) are also 
better represent ed than lower-mass disks. Figure 3 of 
IRobitaille et al.l indicates the presence of "coupled" non- 
uniform sampling between parameters: there are many 
more models with an envelope accretion rate greater than 
10~ 6 M Q yr _1 and that these models predominantly 
have disk masses greater than 10 -3 M Q . 

However, our disk mass estimates are not undermined 
by non-uniform sampling. The lower two panels of 
IRobitaille et all s Fi gure 3 and our Figure [20] indicate 
that simply including far-IR and submm data eliminates 
vast swaths of parameter space and overcomes these sam- 
pling biases. The intrinsic sampling of the Robitaille grid 
is weighted towards massive disks, yet we find numerous 
homologously depleted disks inspite of this bias because 
far-IR/submm non-detections rule out nearly all massive 
disks. Conversely, the submm detections for primordial 
disk s provide strict limits on the inferred disk mass (see 
also IBeckwith et al.|[l990l lAndrews and Williams! l2005ft 
and make the interquartile range of M<ji s fc/M* for these 
objects nearly single valued in many cases (see Figure 
I13[) . On the other hand, uncertainties in deriving other 
disk parameters (e.g. flaring power) have not been fully 
explored here and may be seriously affected by non- 
uniform sampling. 

Our analysis leverages not just on the precision with 
which we can estimate the mass of disks but also our ac- 
curacy. Estimating a mass of emitting dust from submm 
data alone requires assuming a value for the dust opac- 
ity, k (e.s.. IBeckwith et al.lll99C; IBeckwith and Sargentl 
119911: iHenning and Stognienkol Il996l . and later refer- 
ences) . Since the dust opacity is not strictly known (and 
is affected by grain growth), what can be constrained 
from such data is the product of the dust mass and 



the opacity, not simply the dust mass. On the other 
hand, grain growth/dust settling (which affect n) have 
far stronger effects on the near-to- mid IR port i on of 
disk SEDs than in t he submm (e.g. IWood et al.l 120021: 
ID'Alessio et al.ll2006T). For inst ance, "settled disk" mod- 
els from ID'Alessio et all (|2006l ) have 250-1000 [im emis- 
sion equal to or larger than that for disks without grain 
growth/settling (see Figure 13, top two sets of panels in 
ID'Alessio et^l200fih . 

The total disk mass is also uncertain. Since contin- 
uum IR to submm emission comes predominantly from 
dust, a gas-to-dust ratio must be assumed in converting 
a submm-inferred d ust mass to a total disk mass. The 
mod els used here. inlCurrie et al.l (120091) . and other w ork 



(e.g. ICieza et al.H2010t lAndrews and Wnliamsll2005l ) as- 
sume st andard values from k and a solar gas-to-dust ra- 
tio f see IBeckwith and Sargentlll991ft . However, compar- 
isons with disk masses derived from accretion rates (as- 
suming steady-state accretion) indicate that our meth- 
ods may systematically underesti mate true disk masses 
(e.g. lAndrews and Williams! 120071 ). This may be due ei- 
ther to an inaccurate assumed value for k or a substan- 
tially non-solar gas-to-dust ratio, which may occur from 
grain growth and planetesimal formation. Our analysis 
provides evidence that masses for homologously depleted 
transitional disks are systematically lower by a factor of 
10-1000 compared to primordial disks; other work pro- 
vides evidence that gas and small dust deplete on similar 
timescales for disks in general and for homologously de- 
pleted transitional disks in particular (iFedele et aLll2010l : 
Kennedy and Kenv on 2009; Cieza et aLll2010D . However, 
given the uncertainties in determining k and the gas- 
to-dust ratio, our analysis does not provide an abso- 
lute calibration for disk ma sses that can inform mod- 
els for planet formation (e.g . lKenvon and Bromlevll2009l : 
iBromlev and Kenvodl2l)llT . 

Accurately quantifying uncertainties in stellar mass for 
each star is another challeng e. Between 1 Myr and 5 
Myr, the iBaraffe et al.l (|1998h isochrones with a mixing 
length parameter of L p = 1.9 H p ("required to fit the 
Sun") predict that 1.4 M Q stars increase in T e ff by ~ 
125 K from 4663 K to 4786 K , thus changing in spec- 
tral type by ~ 1 subclass (cf. ICurrie et al.l l2010). The 
variation is smaller (~ 100 K) for solar-mass stars and 
inconsequential for very low- mass stars (M* < 0.7 M Q ). 
Adopting the conversion from spectral type to mass from 
IBaraffe et all our stellar masses are then intrinsically un- 
certain by ~ 0.1 M Q or less due to the different mapping 
between spectral type and mass as a function of age. 
The uncertainties in measured spectral types are up to 
1-2 subclasses, considering the range of spectral types 
repo rted in the literature for var ious Taurus members 
(e.g. IWhite and Hillenbrand! 12004 . The age and spec- 
tral type uncertainties listed above indicate that our frac- 
tional disk masses could be uncertain by an additional 
10-20%, and that considering these uncertainties could 
slightly broaden the interquartile range of fractional disk 
masses for some objects. 

Finally, our analysis may miss some sources whose gaps 
cannot be inferred from photometric data alone. Among 
our sample, the SED for DH Tau shows evidence for a 
XF\ oc A -3 decline from 3.6 /im to 8 /j,m, similar to that 
for transitional disks UX Tau and LkCa 15, which have 
large (tens of AU) gaps separating a small optically-thick, 
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~ 1500 K inner disk from an optically-thick outer disk. 
The SED fit to this source from the Robitaille et al grid 
is poor. Since the Robitaille et al grid provides good 
sampling for full disks with a range of masses (and thus 
mid-IR emission) and disks with cleaned inner holes, it is 
plausible that DH Tau's disk has a different morphology 
such as a disk with a developing gap. Since not all Taurus 
members have been targeted by mid-IR spectrographs, 
more disks with small holes / gaps may be id e ntified by 
future observations. Indeed, lEspaillat et al.l (|2011l ) ar- 
gue that IP Tau, which we identify as having a primor- 
dial disk in spite of its weak mid-IR emission, has a ~ 
2 AU gap. Thus, the true size of the transitional disk 
population in Taurus, especially the subset with small 
gaps/inner holes, is plausibly larger than what we find. 

5.4.2. Limits on Our Study As a Probe of Transitional 
Disk Lifetimes/Disk Clearing 

While our analysis provides strong evidence in favor 
of a long transitional disk phase, especially for very 
low-mass stars, several important caveats qualify this 
conclusion. First, many different mechanisms may be 
responsible for explaining the morphologies of transi- 
tiona l disks ([Alexander and Arm itagc 20Q3; ICieza et al.l 
120101 ISicilia-Aguilar et al.1 12010ft . Therefore, a "transi- 
tion timescale/transitional disk lifetime" derived from 
the relative frequency of all transitional disks is only an 
averaged value of the timescales from a number of differ- 
ent mechanisms like gas giant planet formation and UV 
photocvaporation. Some of these mechanisms may oper- 
ate on fast, < 0.1-0.5 Myr timescales but less frequently 
whereas others take > 1 Myr to clear dust in disks and 
operate more frequently. 

The disk clearing time for many candidate mechanisms 
(e.g. gas giant planet formation, UV photoevaporation) 
depends on properties such as disk viscosity, initial disk 
mass and initial angular momentum that are poorly con- 



strained and likely have a large intrinsic dispersion. This 
dispersion in disk properties may produce a large dis- 
persion in the duration of the transition phase similar 
to the observed dispersion in total protoplanetary disk 
lifetimes. While the mean tr ansition time s cale m ay be 
~ 1 Myr as found here, in ICurrie et al.1 (|2009( 1 , and 
ISicilia-Aguilar et "ail (|2009f ). the transitional disk phase 
for individual disks may last much shorter or longer. 

Quantitatively assessing how transitional disks probe 
the clearing of both gas and dust requires diagnostics of 
circumstellar gas sensitive to the bulk gas content. The 
Herschel Space Observatory offers a sensitive probe of 
far-IR line emission to identify cool gas in planet-forming 
regions of the disk. Herschel programs such as GASPS 
will survey many 1-10 Myr-old stars for evidence of cir- 
cumstellar gas and thus provide strong constraints on 
gas dissipation as a function of time. Comparing the gas 
properties of transitional disks with those for primordial 
disks will more definitively determine how transitional 
disks are clearing their gas, complementing studies of 
dust clearing investigated in our work. 
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TABLE 1 
Members Listed in Forbrich and 
Preibisch (2007) And Included in This 
Work 



Name RA (2000) DEC (2000) 
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08 


.3 


IRS7w 


19: 


01 


:55. 


.31 


-36: 


57: 


22. 


.0 


FP-33 


19: 


01 


:55. 


61 


-36: 


56: 


51. 


.1 


FP-34 


19: 


01 


:55. 


76 


-36: 


57: 


27. 


.7 


FP-35 


19: 


01 


:55. 


.85 


-36: 


52: 


01. 


.3 


IRS-7e 


19: 


01 


:56. 


.39 


-36: 


57: 


:28. 


.4 


FP-37 


19: 


01 


:57. 


16 


-37: 


03: 


11 


.9 


FP-38 


19: 


01 


:58. 


.32 


-37: 


00: 


.27. 


.5 


T-CrA 


19: 


01 


:58. 


79 


-36: 


57: 


50. 


.1 


V702 


19: 


02 


:01 


.92 


-37: 


07: 


13 


.0 


B1858 


19: 


02 


:01 


.91 


-36: 


51: 


00 


1 


HBC-679 


19: 


02 


:22. 


.13 


-36: 


.55: 


11. 


.0 


CrA-159 


19: 


02 


:33. 


.07 


-36: 


.58: 


21. 


.1 



NOTE. — Sources with an "FP" prefix followed by 
a number identify Coronet Cluster members uniquely 
identified by Forbrich and Preibisch I 2007)- The 
number following "FP" identifies the row number of 
the source in Table 2 of that paper. 



HD176386B 




TY-CrAabcd 



Fig. 1. — Region of the IRAC [5.8] mosaic showing nebular emission surrounding TY CrA and HD 176386b. The positions of both 
sources are shown as green circles. Nebular emission contaminates both stars and is responsible for positional offsets between the 2.2-4.5 
fim centroids and 5.8—8 fim centroids. 
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TABLE 2 

Photometry for Members Listed in Forbrich and Preibisch (2007) 



Name 


J 


a(J) 


H 


cr(H) 


K s 


<r(K s ) 


[3.6] 


"([3.6]) 


[4.5] 


a([4.5]) 


[5.8] 


ct([5.8]) 


[8] 


"([8]) 


[24] 




S-CrA 


8.19 


0.02 


7.05 


0.024 


6.10 


0.02 


5.060 


0.005 


4.386 


0.005 


3.699 


0.005 


2.969 


0.004 


- 


- 


FP-3 


- 


- 


- 


- 


- 


- 


13.971 


0.014 


13.371 


0.012 


12.092 


0.020 


11.616 


0.029 


- 


- 


FP-6 


- 


- 


- 


- 


- 


- 


- 


- 


- 


- 


12.299 


0.017 


11.458 


0.033 


6.106 


0.0 


FP-8 


- 


- 


- 


- 


- 


- 


- 


- 


- 


- 


- 


- 


12.697 


0.092 


- 


- 


CrA-134 


11.61 


0.02 


9.78 


0.03 


8.77 


0.023 


7.876 


0.002 


7.367 


0.002 


6.862 


0.003 


6.537 


0.002 


- 


- 


CrA-135 


10.81 


0.02 


10.11 


0.02 


9.89 


0.02 


9.599 


0.005 


9.601 


0.004 


9.551 


0.007 


9.611 


0.007 


- 


- 


V709 


8.67 


0.03 


7.97 


0.04 


7.71 


0.02 


7.863 


0.002 


7.717 


0.003 


7.533 


0.004 


7.512 


0.004 


6.794 


0.0 


I HD 


7.60 


- 


7.49 


0.09 


7.14 


0.08 


7.837 


0.003 


7.781 


0.003 


6.243 


0.003 


4.706 


0.003 


- 


- 


176386 


































FP-18 


11.34 


0.03 


10.66 


0.03 


10.37 


0.03 


10.110 


0.005 


10.097 


0.005 


9.979 


0.007 


9.927 


0.010 


9.364 


0.1 


X TY CrA 


7.49 


0.02 


6.97 


0.03 


6.67 


0.02 


6.217 


0.005 


6.085 


0.005 


4.478 


0.004 


- 


- 


- 


- 


abed 


































IRS-2 


13.75 


0.02 


9.74 


0.02 


7.07 


0.02 


5.181 


0.004 


4.281 


0.003 


- 


- 


2.230 


0.003 


- 


- 


HBC 677 


10.47 


0.02 


9.04 


0.02 


8.09 


0.03 


7.725 


0.003 


7.201 


0.003 


6.528 


0.004 


5.859 


0.004 


- 


- 


FP-23 


16.92 


0.18 


15.23 


0.11 


14.77 


0.10 


- 


- 


- 


- 


- 


- 


- 


- 


- 


- 


IRS-5ab 


16.9 


- 


13.73 


0.01 


10.59 


0.01 


7.456 


0.002 


6.293 


0.001 


- 


- 


3.937 


0.003 


- 


- 


FP-25 


15.56 


- 


13.14 


- 


12.49 


0.15 


- 


- 


9.115 


0.003 


- 


- 


7.017 


0.003 


- 


- 


IRS-6A 


16.39 


0.20 


12.14 


- 


10.31 


- 


8.480 


0.003 


7.687 


0.004 


- 


- 


6.286 


0.003 


- 


- 


V710 


15.96 


0.10 


11.13 


0.03 


7.62 


0.04 


- 


- 


3.659 


0.003 


- 


- 


- 


- 


- 


- 


IRS-8 


14.92 


0.04 


12.69 


0.03 


11.37 


0.02 


9.911 


0.005 


9.358 


0.005 


8.744 


0.006 


7.946 


0.005 


- 


- 


IRS-9 


- 


- 


- 


- 


- 


- 


8.031 


0.003 


7.313 


0.003 


- 


- 


- 


- 


- 


- 


R CrA 


6.952 


0.021 


4.955 


0.024 


2.910 


0.318 


- 


- 


- 


- 


- 


- 


- 


- 


- 


- 


IRS-7w 


12.55 


- 


13.25 


0.32 


10.76 


- 


8.635 


0.004 


7.073 


0.003 


- 


- 






- 


- 


FP-33 














10.181 


0.005 


9.025 


0.004 






6.860 


0.004 






FP-34 


- 






- 


- 


- 


9.730 


0.004 


8.293 


0.004 














r r-oO 


















13.387 


0.016 


ll.DDU 


U.Ulo 




U.UzU 


I .ooZ 


n n 


IRS-7e 














9.417 


0.004 


















FP-37 


15.58 


0.07 


14.02 


0.04 


13.19 


0.04 


12.320 


0.007 


12.264 


0.006 


11.925 


0.089 


12.086 


0.037 






FP-38 


16.961 




15.275 


0.139 


13.704 


0.060 


12.426 


0.007 


12.102 


0.008 


11.508 


0.017 


12.008 


0.050 






T-CrA 


8.93 


0.03 


7.70 


0.04 


6.60 


0.02 






4.819 


0.005 


4.217 


0.005 










V702 


8.90 


0.02 


8.48 


0.04 


8.35 


0.03 


8.444 


0.004 


8.329 


0.003 


8.282 


0.0048 


8.280 


0.003 


8.122 


0.0 


B185839.6 














14.099 


0.019 


14.766 


0.054 










8.617 


0.0 


-3658 


































HBC-679 


10.33 


0.03 


9.50 


0.04 


9.23 


0.03 


9.031 


0.003 


8.967 


0.004 


8.901 


0.005 


8.868 


0.006 


8.258 


0.0 


CrA-159 


10.59 


0.03 


9.26 


0.02 


8.45 


0.020 


7.565 


0.002 


7.016 


0.003 


6.372 


0.003 


5.887 


0.004 


3.704 


0.0 


Note. — Notes 


- 1) Inspc 


ction of the IRAC n 


rosaic shoi 


vs that TY 


CrAabcd z 


nd HD 176386B arc cor 


taminatcd by nebular e 




5.8 fim and 


3 fm. 2). 









Fig. 2. — (Left panel) I RAC color-color diagram of Coronet Cluster stars with new photometry (black circles) and those analyzed in 
ISicilia- Aguilar et al. (2008) (grey circles). (Right panel) K s -[5.8] vs. K s -[24] color-c olor diagram for C o ronet Cluster stars. In each panel, 
we overplot a r edden ing vector of Ay = 10, assuming the extinction relations of Indcbctouw ct al. (2005) for 2MAS S and IRAC and 



Flaherty ct al. (2007) for MIPS. In the left panel, the box corresponds to the range of colors for most stars in Taurus (Hartmann ct al. 
2005); the dotted box in the right-han d panel corresponds to the range of colors for most homologously depleted transitional disks in NGC 
2362 identified by Curric ct al. (2009). The triangles in the right panel corre spond to (from bott om to top) the lower-quartile Taurus SED, 
median Taurus SED, and upper- quartile Taurus SED for K5-M2 stars from lFurlan et al.l (120061, . 
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TABLE 3 

Provisional Evolutionary States for Coronet Cluster Disks with New Spitzer Photometry 

Based on 1-24 jim Data 



ID 


ST 


ST Ref 


A v (best) 


T mid-IR 


Inner Hole? (R, 


sub) Disk Class 




vjO/ xvO 


i 


9 


I111CK 


ii 


p n 
i \j 


CrA-159 


M2 2 


est. 


3 


thick 


n 


PD 


IRS-2 


K2 


2 




thick 


n 


I >c 

la 


CrA-134 


K3 


2 


11 


thick 


n 


PD 


CrA-135 


M4 


2 


0.5 


— . 




Star 


v /uy 




^7 
( 


0.5 


thm/— 




otar/ 1 D(riD J / : 




K7 


2 


0.75 




n 




FP-18 


M3 


est. 


0.4 






Star 


TY-CrA 


B8/K2 


. > . i 




? 




7 


HBC-677 


M2 


2 


3.9 


thick 


n 


PD 


IRS-5ab 


K5 


2,5 




thick 


n 


PS 


FP-25 


M? 


7 




thick 


n 


PS 


IRS-6 


Ml 


2 


10.5 


thick 


n 


PD 


V710 


K5-M0 


5 




thick 


n 


PS 


IRS-7 


? 






thick 


n 


PS 


IRS-8 


M2 


est. 


13.5 


thick 


n 


PD 


R-CrA 


A5 


7 


? 


thick 


n 


PD 


T-CrA 


FO 


6 


0.3 


thick 


n 


PD 


V702 


G5 


7 


0.25 






Star 


HBC-679 


K2 


7 


1.5 






Star 


Note. - 


The disk states are 


identified 


as follows: 


PS=Protostar, 


PD = Primordial Disk, 



and TD = Transitional Disk, Star = Stellar Photosphere (no circumstellar material). The 
transitional disks are further divided into those with inner holes (IH) and those that are homol- 
ogously depleted (HD). The disk states here are considered to be provisional since these Coronet 
Cl uster stars lack s ensitiv e fa r-IR/submm data. The spe ctral type references are the following: 
1) iCarmona et all (|2007l). 2) iMever and Wilkind (l2009f> 3) [Marraco and Rvdgrenl (119811), 4) 
ICasevet all l|1998l ). 5) iNisini et alj l|2005D . 6) iPattenl (|1998l ). and 7) iForbrich and Preibischl 
(2007). "Est" means that the spectra l type was esti mated by modeling the optical to near-IR 
SED using the iRobitaille et all (|2006T ) grid and the iCurrie et al.l (|2010l ) effective temperature 
scale. Note (1) - The fit to the photosphere for this star is generally poor; the best-fit Av 
corresponds to a K0 stell ar photospher e model. However, the exact value has no bearing on our 
classification. Note (2) - Patten ( 1998) list the spectral type for this star as M3-M5. However, 
we canno t reproduce the ob served optical to near-IR SED with a M3-M5 photosphere using 
either the lCurrie et al.l (|201lD Te scale and dwarf colors or with the lRobitaille et al.l (|2006l ) grid: 
all good-fitting models (x 2 — Xbest < 3) predict that its photosphere is hotter than ~ 380 K and 
thus i ts spectral type is earlier than M0-M2. Since other stars classified as M3-M5 by|Pattcn 
l|1998D are reclassified as M0-M2 by IMever and W ilking (2009) based on higher signal-to-noise 
data, we list this star as an M2 star. 
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TABLE 4 

Provisional Evolutionary States for Coronet Cluster Disks Studied in Sicilia-A guilar et al.I (120081 1 Based on 

1-70 fim Data 



ID 


ST 


ST Ref. 


A v (best) 


7~mid — I R 


Inner Hole?(R su f,) 


Disk Class 
(Si08) 


Disk Class 
(Er09) 


Disk Class 
(This Work) 


/"!„ A rtr\ C 

OrA-2U5 


~\ ft A TV If C 

1V14— Mo 


8,9 


0.25 


. 

trim 


y,2000 


rnn /T"LJ\ 

1 D(lri) 


PD 


1 D(LH) 


/~i„ A A on 

OrA-4oz 


Mo— M / 


8,9 


0.5 


thick 


n 


ID 






A /ICQ 

OrA-4oo 


M4.5 


10 


1 . 5 






otar 




btar 


A /ICC 

OrA-4oo 


TV K C TV J'7 C 


8,10 





thick 


n 






JrU 


/~V A /ifif; 
v_ 1V\-40D 


1\1Z 


1 n 

J.U 


D.u 


thick 


n 


1 D (111 J 


1 LJ 


pn 
1 L> 




>* 1V14.0 


S Q 


1 


thick 


n 


pn 








.> 1V10 


Q oof 

y , esr 


1 

1 


L 111 11 


y,ouu 






_L U (in j 


CrA-4110 


M5 


9 


O.i 


thick 


n 


PD 




PD 


CrA-4111 


M4.5 


9 





thin 


y,70 


TD(IH) 


TD(IH) 


TD(IH) 


G-l 


MO 


9, est. 


3 


thin 


2 


PD 




TD(HD) 


G-14 


M4.5 


10 


1.9 


thin 


n 


TD(IH) 


PD 


TD(HD) 


G-30 


M3.5 


10 


0.09 


thin 


n 


TD(IH)/DD 




TD(HD) 


G-32 


> M5 


10,est. 


8-15 


thick 


n 


PD 




PD 


G-49 


M4 


10 


0.07 






Star 




Star 


G-65 


M1-M2 


10 


14-15 


thick 


n 


TD(IH) 


PD 


PD 


G-85 


M0.5 


9, est. 


17 


thick 


n 


PD 




PD 


G-87 


Ml. 5 


10 


11 


thin 


1 


TD(IH) 


PD 


TD(HD) 


G-94 


M3.5 


10 


0.6 


thin 


n 


TD(IH)/DD 




TD(HD)/DD 


G-95 


Ml 


10 


5 


thin 


y,20 


TD(IH)/DD 




TD(IH)/DD 


G-102 


M5 


10 


0.7 


thin 


n 


TD(IH)/DD 




TD(HD)/DD 



NOT E. — References for spectral types arc the following: 8) |L6p6z-Marti et aT7"| (20051 ) , 9) Sicilia-Aguilar et al. 2 010, in prep., 10)|SiciUa-Aguilar et "aTl 
120081) . "Est" mea ns that the spectral type was estimated by modeling the optical to near-IR SED using the IKobitaUle et alj <20uTl) grid and the 
ICurrie et alj 120101 ) effective temperature scale. 



TABLE 5 ^ f 

IRS Flux Density Measurements Used in Modeling Coronet Cluster Stars from [Sicilia-Aguilar et al.I (120081) 



ID 


Wavelengths Sampled (|im) 


Flux Density (mjy) 


CrA-205 
G-14 
G-65 
G-87 


9,13,17,21,27 

6,7,9, 10, 11, 12 
9, 13, 17, 21, 27 


1.03 ± 0.06. 0.70 ± 0.08. 0.91 ± 0.30. 1.10 ± 0.30, 2.08 ± 0.26 
4.95 ± 0.20, 4.70 ± 0.15, 5.53 ± 0.45, 4.75 ± 0.34, 3.62 ± 0.49 

53.06 ± 12.80, 48.30 ± 30.30, 49.96 ± 13.27, 47.9 ± 12.30, 48.10 ± 23.98, 43.1 ± 32.70 
10.02 ± 0.31, 12.56 ± 0.29, 15.55 ± 4.84, 19.04 ± 4.88, 25.37 ± 6.19 



NOTE. — The wavelength sampling for G-65 is different than the others because G-65 lacks data longwards of ~ 15 fxm 



TABLE 6 

Properties of Selected Taurus Sources 



Name 


ST 


M* 


A v 


Ay range 


K s -[8] 


K s -[24] 


T mid-IR 


Inner 






(M ) 












Hole? 


TTY Tan 


K r i 
xvo 


i 


o 


u. / o 


fi 9c; 1 OK 


1 C\R 


5 42 


thin 


y 


VS07 Tan 


K5 


i 





1.4 


0.9-1.9 


1.46 


3.91 






T.kflfl 1 ^ 

TjIY V_y a lij 


K5 


i 


o 


1 2 


0.7-1.7 


1.58 


5.06 


thin 


y 


FY Tau 


K5 


l 





5 


4—6 


1.68 


3.93 


thin 




JH 112 


K6 





95 


3.4 


2 72-4 08 


1.97 


5.34 


thick 




VS^fi Tan 

V OOU lull 


K7 


o 


85 


1.05 




1.76 


4.77 






DN Tau 


MQ 


o 


75 


0.7 


0.2-1.2 


1.99 


5.00 


thin 


n 


VY Tau 


MO 


o 


75 


1.1 


0.6-1.6 


1.53 


4.35 


thin 


n 


CIO Tan 


MO 


o 


75 


2 


1.5—2.5 


1.85 


4.93 




11 


IP Tau 


MO 


o 


75 


u.o 


1—1 1 


l.UO 


4 74 


thick 




HO Tan 


MO. 5 


o 


65 


1 3 


0.8—1.8 


1.90 


4.79 


thick 


J| 


F)H Tan 


Ml 


o 


(-; 


1.9 


1.4—2.4 


1.28 


4.79 




? 


FX Tau 


Ml. 5 


o 


57 


2 


1.5—2.5 


1.90 


4.85 


thick 




CY Tau 


Ml. 5 


o 


57 


1.7 


1.2—2.2 


1.89 


4.11 


thick 






Ml. 75 


o 


57 


7 9 


2.4—3.6 


1.49 


5.81 


thin 


y 


-1-9997907 




















TH-223 


M2 





55 


2.2 


1.6—2.6 


1.72 


4.27 


thin 




HH Tan 


M2 





55 




0.5—1.5 


1.79 


4.51 


thick 




ZZ Tau 


M3 





1 


1.2 


0.7-1.7 


1.39 


3.91 


thin 


11 


XEST13-10 


M3 





4 


6 


4.8-7.2 


1.52 


4.40 


thin 


n 


CIDA 8 


M3.5 





3 


2.25 


1.75-2.75 


1.49 


4.30 


thin 


11 


04231822 


M3.5 





3 


8 


6.4-9.6 


1.53 


5.11 


thin 


11 


+2641156 




















04153916 


M3.75 





275 


2 


1.5-2.5 


1.62 


4.93 


thin 


n 


+2818586 




















FP Tau 


Ml 





25 


0.5 


0-1 


1.52 


4.53 


thin 


11 


ITG-15 


M5 





175 


3.25 


2.6-3.9 


1.62 


4.79 


thin 


11 


04202555 


M5.25 





15 


2.1 


1.6-2.6 


1.69 


5.37 


thin 


? 



+2700355 



NOT E. — The K s -[8] and K s -[24] colors listed are deredden e d colo rs as determined by I Lull man et al.l 
<2010|) . The stel lar masses are c ompu ted using the |Currie et al.| <2011|) effective temp erature scale for K5- 
M2 stars and the lLuhman et aU 12003) scale for later stars and the lBaraffe et al.l 1 19 9 SI ) isochro nes. A v was 
determ ined individually by fitting the optical-to-ncar-IFt SED for stel lar photospheres usi ng the|CurHe_e^aL 
I2(mj) intrinsic colors for a given effective temperature and using the lCurrie et all 1201 U) and |LuhirTaTTTFTr 
120*03) effective temperature scales. The range in Ay assumes a nominal 20% extinction uncertainty with a 
minimum uncertainty of 0.5 mags. 
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TABLE 7 

Disk Evolutionary States for Selected Taurus Sources 



Name 


IVl rlo o t> 


M di3k 
M„ 




MdiBk 

(TOR") 


f(< 1,3 xl0~ 3 ) 


Submm. 
Ref. 


a,is ft 
V 0, submm) 


M t 

( qi ibm m 1 

1 5UU111111 j 


Disk Class 


UX Tau 


0.02 


0.02 




[0.017,0.022] 


0,0 


1 


5.1xl0" 3 


5.1xl0" 3 


TD(IH) 


V807 Tau 


3xl0" 3 


3xl0" 3 




[3xl0- 3 ,4xl0" 3 ] 


0,0.5 


1 


10~ 3 


10" 3 


PD/TD(HD) 


LkCa 15 


0.01 


0.01 




[0.01,0.025] 


0,0 


1 


0.05 


0.05 


TD(IH) 


FY Tau 


io- 4 


10" 4 




[6xl0" 6 ,lxl0^ 3 ] 


0.76,0.82 


1 


7xHT 4 


7xl0~ 4 


TD(HD) 


JH 112 


8xl0~ 3 


9.4x10" 


3 


[4xl0" 3 ,0.01] 


0.06,0.12 


1 


10" 3 


1.05xl0~ 3 


PD 


V836 Tau 


0.01 


0.012 




[0.012,0.014] 


0,0 


1 


0.01 


0.012 


PD 


HM Ton 


U.Ul 


U.Ulo 




rn ni q n nQl 
[U.Ulo, U.UoJ 


n n 
u,u 


"1 1 
1.1 


OA 1U 


,\ w i n — 2 
X 1U 




VY Tau 


2xl0~ 4 


2.7x10" 


4 


[3.4xl0~ 5 ,9xl0" 4 ] 


0.76,0.91 


1 


< 5xl0" 4 


<6.7xl0~ 4 


TD(HD) 


GO Tau 


6xl0~ 3 


8xl0~ 3 




[5xl0" 3 ,0.023] 


0,0 


1 


0.07 


0.09 


PD 


IP Tau 


8xl0~ 3 


0.011 




[7xl0- 3 ,0.016] 


0,0.05 


1 


3xl0" 3 


4xl0~ 3 


PD 


HO Tau 


8xl0~ 3 


0.012 




[0.013,0.013] 


0,0 


1 


2xl0" 3 


3.1xl0" 3 


PD 


DH Tau 


5xl0" 3 


0.01 




[0.01,0.018] 


0,0.27 


1 


3xl0" 3 


5xl0" 3 


PD/TD(IH) 


FX Tau 


4xl0~ 3 


7xl0" 3 




[5xl0" 3 ,0.013] 


0.09,0.22 


1 


9xl0" 4 


1.5xl0" 3 


PD 


CY Tau 


8xl0~ 3 


0.014 




[0.014,0.025] 


0,0 


1,4 


6xl0" 3 


0.015 


PD 


04333905 


4xl0~ 4 


7xl0~ 4 




[l.7xl0" 4 ,2.7xl0- 3 ] 


0.63,0.68 


_ 


_ 


_ 


TD(IH) 


+2227207 


















JH 223 


7.9xl0 -5 


1.4x10" 


4 


[8.5xl0- 5 ,lxl0" 3 ] 


0.75,0.82 


1 


<3xl0~ 4 


<5.5xl0- 4 


TD(HD) 


GH Tau 


4.1xl0" 3 


7.5x10" 


3 


[4xl0- 3 ,0.01] 


0.03,0.16 


1 


7xl0" 4 


1.3xl0~ 3 


PD 


ZZ Tau 


lxl0~ 4 


2.5x10" 


4 


[2xl0~ 5 ,8xl0~ 4 ] 


0.77,0.86 


1,3 


<4xl0~ 4 


<10" 3 


TD(HD) 


XEST13-10 


4xl0~ 4 


lxlO" 3 




[2.5xl0~ 4 ,3xl0" 3 ] 


0.5,0.72 








TD(HD) 


CIDA 8 


8xl0" 3 


0.03 




[0.03,0.04] 


0,0 


1 


10" 3 


3.3xl0" 3 


PD 


04231822 


5xl0" 4 


1.7x10" 


3 


[3xl0" 4 ,5xl0^ 3 ] 


0.43,0.66 








PD/TD(HD) 


+2641156 




















04153916 


4xl0~ 4 


1.5x10" 


3 


[4xl0" 4 ,6xl0- 3 ] 


0.39,0.59 








PD/TD(HD) 


+2818586 




















FP Tau 


2.5xl0" 4 


lxlO" 3 




[4xl0" 4 ,4xl0" 3 ] 


0.53,0.68 


2 


<1.3xl0" 4 


<5.2xl0" 4 


TD(HD) 


ITG-15 


3xl0~ 4 


1.7x10" 


3 


[3.6xl0- 4 ,9xl0" 3 ] 


0.49,0.65 








PD/TD(HD) 


04202555 


4xl0~ 4 


2.7x10" 


3 


[Ixl0- 3 ,3xl0^ 3 ] 


0.07,0.53 








PD/TD(IH) 



+2700355 



Not e. — IQR refers to the interquartile range of fractional di sk masses from best-fi t models. The identification "LulO" refers to lLuhman et all 
where the classification listed is based on applyin g the [Luhmari^^l. color criteria. The references for submillimctcr data and 

(su b)millimeter-derived disk ma sses are the following: (1) lAndrews and Williams] lhq05|). f2HJewitt et al.1 lll994f ). (3) IJensen et"aTl (I1994D . and 
(4) IBeckwi th a^id_Sar^entjJ|1991 [) ■ T he fractional disk masses are computed using thefCurrie et^dTl201lT) ~effective temperature scale for K5— M2 
stars and the u^ulrmgmet'all 120031) scale for later stars and the IBaraffe et all fll998l) isochrones. Sources with disk evolutionary states listed 
as PD/TD(IH) or PD/TD(HD) have uncertain states: cither SED modeling is inconclusive as to whether the disks likely have inner holes 
(PD/TD(IH)) or the classifica t ion d epends on whether the full range or interquartile range of disk masses is used to identify primordial disks 
(PD/TDfHD)). ILuhman et al.1 ( I201CD label all of those disks as primordial disks. 



TABLE 8 

Optically- Thick Primordial Disks Around K5-M6 Stars in Taur us Listed As Having Low Masses In 

IAndrews and Williams (20051) 



Name 


ST 


M* 


M dtsk (M Q )(AW05) 


M dlsfc /M, (AW05) 1 




(Mq, this work) 


M disk 


/M 


HN Tau 


K5 


1.0 


8xl0 -4 


8X10"" 4 


6x10" 


3 


6x10 


-3 


V955 Tau 


K7 


0.85 


5xl0" 4 


5.9xl0" 4 


4x10" 


3 


5x10- 


-3 


DF Tau 


M2 


0.55 


4xl0" 4 


7.3xl0" 4 


5x10" 


3 


9x10" 


-3 


CoKu Tau/3 


Ml 


0.6 


<4xl0" 4 


<6.7xl0" 4 


5x10" 


3 


8x10" 


-3 


DP Tau 


Ml 


0.6 


<5xl0" 4 


<8.3xl0" 4 


2x10" 


3 


3x10- 


-3 


CIDA-3/V410 X-Ray 1 


M2 


0.55 


<4xl0" 4 


<7.3xl0" 4 


2x10" 


5 


4x10" 


-5 


CZ Tau 


M3 


0.4 


<4xl0" 4 


<10" 3 


2x10" 


5 


5x10" 


-5 



(this work) 



Note. — 



Notes (1) — Fractional disk mass is derived using our estimates for stellar mass. 
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TABLE 9 

Provisional Disk Evolutionary States for Selected Late-Type IC 348 Members based on 1-24 (im Data 



















Disk Class 




ID 


ST 


A„(Best) 


K s -[8] 


K s -[24] 


'mid — I R 


Inner Hole? 


CK09 


LulO 


This Work 


21 


K0 


5.4 


1.03 


5.43 


thin 


y 


TD(IH) 


TD 


TD(IH) 


26 


K7 


8 


2.03 


4.90 


thick 


n 


PD 


PD 


PD 


32 


K7 


4.6 


1.95 


4.81 


thick 


n 


PD 


PD 


PD 


40 


K8 


4 


1.58 


5.68 


thin 


? 


PD 


PD 


TD(HD)/TD(IH) 


55 


M0. 5 


11 


1.93 


4.86 


thick 


n 


PD 


PD 


PD 


58 


Ml. 25 


3.7 


1.81 


5.81 


thin 


? 


PD 


PD 


TD(HD)/TD(IH)? 


67 


M0. 75 


1.2 


1.14 


6.12 


thin 


y 


TD(IH) 


TD 


TD(IH) 


68 


M3.5 


2.25 


1.43 


4.14 


thin 


? 


PD 


PD 


TD(HD)/TD(IH)? 


72 


M2.5 


1.7 


0.75 


5.65 


thin 


y 


TD(IH) 


TD 


TD(IH) 


76 


M3.75 


3.0 


1.48 


4.41 


thin 


n 


PD 


TD 


TD(HD) 


97 


M2.25 


5.5 


1.30 


4.85 


thin 


y 


PD 


PD 


TD(IH) 


110 


M2 


5.0 


1.62 


5.90 


thin 


y 


PD 


PD 


TD(IH) 


128 


A12 


2.75 


1.47 


3.54 


thin 


n 


PD 


PD 


TD(HD) 


129 


M2 


2.0 


1.88 


4.32 


thin 


n 


PD 


PD 


TD(HD) 


135 


M4.5 


0.9 


1.41 


3.48 


thin 


n 


PD 


TD 


TD(HD) 


"1 ^ 
1 Oo 


M4.75 


2.85 






LI11CK 


n 


PD 


PD 


PD 


194 


M4.75 


3.1 


1.46 


4.88 


thin 


n 


PD 


PD 


TD(HD) 


213 


M4.75 


1.4 


1.48 


4.04 


thin 


n 


PD 


PD 


TD(HD) 


214 


M4.75 


1.1 


1.48 


4.69 


thin 


n 


PD 


PD 


TD(HD) 


301 


M4.75 


6.2 


1.24 


4.93 


thin 


y 


TD 


TD 


TD(IH) 


308 


M4 


11.8 


1.31 


4.92 


thin 


y 


PD 


TD 


TD(IH) 


1928 


M5.5 


1.4 


1.74 


3.68 


thick 


n 


PD 


PD 


PD 


8078 


M0. 5 


5.8 


1.82 


5.22 


thick 


n 


PD 


PD 


PD 


9024 


M0 


3.75 


1.94 


4.72 


thick 


n 


PD 


PD 


PD 


10352 


Ml 


3.5 


1.88 


5.24 


thick 


n 


PD 


PD 


PD 


Note. — 


- Ay (Best) corn 


:sponds to the 


extinction 


that produ 


ccs the best 


fit to the observed 


SED. The ra 


,nge in A V used 


for SED modeling with 


the Robita 


ills models follow 


s the method d 


escribed fc 


»r Taurus st 


ars. The idci 


itif'ication " LulO" f 


iters to Luhn 


lan et al.| 12010). 


where the classification 


listed is b£ 


iscd on applying the Luhman ct 


al.| ISOiCh 


color criteri 


a. The classi 


fication listed for CK09, Currie 


and Kenyan 12009). is based on applying 



their color-criteria. The disk states listed here are considered to be "provisional" since IC 348 lacks sensitive far-IR and submm data. 



TABLE 10 

Provisional Disk Evolutionary States for Late-Type NGC 2362 Members Based on 1-24 fim Data 



ID 


ST 


A„(Best) 


K s -[8] 


K s -[24] 


T mid-IR 


Inner Hole? (R su() ,AU) 


CL09 


Disk Class 
LulO 


This Work 


3 


M2 1 


0.6 


1.45 


5.89 


thin 


y,4.10 


TD(IH) 


PD 


TD(IH) 


36 


K5 1 


0.5 


1.70 


4.51 


thin 


n 


TD(HD) 


I'D 


TD(HD) 


ll 2 


Kl 


0.3 


0.19 


1.34 


thin 


n 


TD(HD) 




TD(HD) 


63 2 


M0 1 


0.7 


0.52 


2.76 


thin 


n 


TD(HD) 




TD(HD) 


85 


M2 1 


1.5 


1.55 


4.58 


thin 


n 


TD(HD) 


PD 


TD(HD) 


111 


K7 


0.75 


2.5 


5.34 


thick 


n 


PD 


PD 


PD 


139 


M2 


0.4 


3.48 


6.82 


thick 


n 


PD 


PD 


PD 


168 


K3 


0.5 





3.48 


thin 


y,5.20 


TD(IH) 


TD 


TD(IH) 


177 


M0 1 


0.8 


1.82 


5.33 


thin 


? 3 


TD(HD) 


PD 


PD /TD(HD) /TD(IH)? 


187 


M0.5 


1.5 


2.64 


4.81 


thick 


n 


PD 


PD 


PD 


194 


K5 


1.2 


1.06 


5.46 


thin 


y,7.57 


TD(IH) 


TD 


TD(IH) 


196 


K5 1 


0.1 


0.22 


2.35 


thin 


y,3.63 


TD(IH) 


TD 


TD(IH) 


202 


M2.5 


0.9 


2.53 


5.92 


thick 


n 


PD 


PD 


PD 


204 


M2.5 


2 


1.46 


4.69 


thin 


n 


TD(HD) 


PD 


TD(HD) 


213 


M3 1 


0.8 


1.28 


3.84 


thin 


n 


TD(HD) 


TD 


TD(HD) 


214 


K7 1 


0.35 


1.5 


4.17 


thin 


n 


TD(HD) 


PD 


TD(HD) 


219 


M0 


0.8 


1.65 


3.87 


thin 


n 


TD(HD) 


PD 


TD(HD) 


229 


Ml. 5 


0.95 


1.38 


3.92 


thin 


n 


TD(HD) 


PD 


TD(HD) 


251 


M2 1 


1.25 


1.7 


4.68 


thin 


n 


TD(HD) 


PD 


TD(HD) 


267 


M2 1 


0.85 


1.33 


4.19 


thin 


n 


TD(HD) 


PD 


TD(HD) 


295 


M2 1 


0.4 


0.34 


4.20 


thin 


y,1.29 


TD(IH) 


TD 


TD(IH) 



NOTE. — The identification "LulO" refers to ILuhman et all 120101), where the clas sification listed is based on applying the ILuhtnan et all <201CI) color criteria. The 
identification CL09 refers to the disk evolutionary states listed in |Gurrie et al.| 1 20091) ■ The disk states liste d from this work are considered to be provisional since NGC 
2362 lacks sensitive far-IR/submm data. To compare our results with those presented in |Currie et al.| 120091 ), we list the inner hole sizes in AU. (1) Denotes approximate 
spectral types assigned based on fitting the optical and near-IR. data. Likewise, Ay (Best) is corresponds to the extinction that produces the best fit to the observed 
SED. (2) The detections of these sources were confirmed by a rcrcduction of the MIPS-24 data and PRF-fitting photometry. (3) The best-fit models are evenly divided 

between those that have an AU-scale inner hole and those that don't. For members with disks that have an inner hole (marked "y"), we add the inner hole size in AU as 
determined in |Currie et al.| f2009T ) . 
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TABLE 11 

Disk Evolutionary States for Late- Type r) Cha Sources 



Name 


ST 1 


My, 




K s -[24] 


' mid — 1 ti 


Inner 
Hole? 






M d „ k 


Disk Class 
f(< l,3xl0 -3 ) 


Si09 


LulO 


This Work 


RECX-l 


K7+M0,K6 


0.95 


0.23 


0.38 


- 


n 


- 




- 


- 


star 


star 


star 


RECX-3 


M3.M3.25 


0.35 


0.26 


0.92 


thin 


y 


3.5x10- 


8 


lxlO- 7 


1,1 


TD(IH) 


TD 


TD(IH) 


RECX-4 


M1.75,M1.3 


0.6 


0.30 


0.89 


thin 


y 


1.8x10" 


7 


3xl0- 7 


1,1 


TD(IH) 


TD 


TD(IH) 


RECX-5 


M3.8.M4 


0.25 


0.97 


4.81 


thin 


y 


7xl0~ 4 




2.8xl0" 3 


0.14,0.73 


TD(IH) 


TD 


TD(IH) 


RECX-6 


M3,M3 


0.375 


0.25 


0.41 


- 


n 






- 




star 


star 


star 


RECX-7 


K6.9+M1,K6 


0.95 


0.17 


0.26 




n 










star 


star 


star 


RECX-9 


M4.4+M4.7, 


0.225 


1.37 


3.96 


thin 


n 


1.4x10" 


4 


6xl0- 4 


0.53,0.75 


TD(HD) 


PD 


TD(HD) 




M4.5 






















RECX-10 


M0.3.M1 


0.6 


0.24 


0.28 




n 










star 


star 


star 


RECX-11 


K6.5,K5.5 


1 


1.69 


3.98 


thin 


n 


3.3x10- 


3 


3.3xl0- 3 


0.13,0.50 


PD 


PD 


PD 


RECX-12 


M3.2,M3.25 


0.35 


0.34 


0.46 




n 










star 


star 


star 


J0841 


M5.5,M4.75 


0.2 


1.50 


3.91 


thin 


n 


3xl0" 5 




2xl0- 4 


1,1 


TD(HD) 


PD 


TD(HD) 


J0843 


M5,M3.25 


0.35 


2.92 


5.91 


thick 


n 


lxlO" 3 




3xl0- 3 


0.31,0.50 


PD 


PD 


PD 


J0844 


M5.3, M5.75 


0.1 


1.86 


4.51 


thin? 


n 










PD 


? 


PD/TD(HD) 2 


J0838 


M5, M5.25 


0.15 


0.53 






n 










star 


star 


star 


J0836 


M5.3, M5.5 


0.12 


0.58 






n 










star 


star 


star 



scriptsize 

Note. — The identification "LulO" refers to Luhman ct al. (2010) , w here the classification listed is based on applying the Luhinan ct al. (2010) 
color criteria. The identi fication "Si 9" ref ers to Sicilia-Aguilar ct al. (200j| )'s disk classification t r anslat ed into our nomenclature. 1) The first entry for 
spectral type comes from Lyo ct al. (2004) and the second one comes from Luhman and Stccghs (2004). We compute the stellar masses and fractional 
disk masses using the Luhman and Stecghs (2004) spectral types and the stellar masses for a given effective temperature at 6 Myr from the Baraffc ct al. 
(1998) isochrones. Our results do not leverage on whose spectral types we adopt. 2) For our statistics, we classify this disk as a primordial disk since 
it la cks far-IR data t o needed to derive a meani ngful disk mass estimate. We omit a column for optical extinction since it is negligible for this cluster 
fe.e. lLvo et~aT]|200l ; ILuhman and Steeghsl 120041) . 



TABLE 12 

Frequency of Transitional Disks With Time (Corrected for 
Contamination) 



Cluster 


Age (Myr) 


Spectral Type Range 


f (TD); ou>eT . 


f (TD)^pp er 


Taurus 


1-2 


K5-M6 


0.18 (25/143) 


0.22 (32/145) 






K5-M2 


0.17 (13/75) 


0.21 (16/77) 






M2.5-M6 


0.18 (12/68) 


0.24 (16/68) 


Coronet 


1-3 


K5-M6 


0.21 (4/19) 


0.36 (8/22) 






K5-M2 


0.12 (1/9) 


0.22 (2/9) 






M2.5-M6 


0.30 (3/10) 


0.50 (6/12) 


IC 348 


2.5 


K5-M6 


0.32 (40/124) 


0.40 (49/124) 






K5-M2 


0.18 (6/33) 


0.24 (8/33) 






M2.5-M6 


0.42 (34/81) 


0.51 (41/81) 


NGC 2362 


5 


K5-M3 


0.58 (11/19) 


0.63 (12/19) 


r) Cha 


6-8 


K5-M6 


0.63 (5/8) 





Note. The transitional disk frequency is denned as 

f(TD)/(f(TD)+f(PD)). The range in transitional disk frequencies (f(TD)) 
for Taurus accounts for disks with uncertain/indeterminable states. The 
range in frequencies for the Coronet, IC 348, and NGC 2362 also account 
for the range of contamination estimates based on SED modeling of Taurus 
data. 
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Fig. 3. — Spectral energy distributions of Coronet Cluster stars with new Spitzer photometry. In most cases, the SEDs can be modeled as 
a stellar photosphere (gree n) pl us excess emission from a disk. The stellar photospheres plotted come from the MARCS stellar atmosphere 
models (Gustafsson ct al. 2008, and references therein). Sources not shown lack near-to-mid IR data required to constrain their SEDs. 
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Fig. 4. — Compari s ons b etween various empirical and theoretical SEDs and our fiducial flat disk models used for analy s is- the 
Kcnyon and Hartmann (1987) flat, reprocessing disk and our non-flaring, optically thick reprocessing disk f rom the Whitney c t alj l|2003a| ) 
Monte Carlo radiative transfer code. The models shown are for MO, T e = 3850 K stars. In both panels, the Kcnyon and Hartmann (1987) 
disk is shown as a dashed-three dots— diamonds locus, the [Whitney ct al. model is shown as a grey dotted l ine, and the s tellar photo spher e 
is shown as a solid line connected by squares, (left) Comparisons with the lower-quartile Taurus SED from Furlan ct al. (dash/dots 2006). 
(right) Comparisons with a standard flared disk model (grey dashed line) and a flared b ut " settled" / shadowed d isk model (black dashed 
line). The Whitney radiative transfer models deviate from the stellar photosphere and Kcnyon and Hartmann flat disk in the near-IR 
because it includes strong, hot, optically-thick emission from the puffed up inner disk regions. 
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Settled, Optically Thick Disk 



Flattened, Optically Thick Disk 




Ks-[24] 




Settled, Optically Thick Disk 



Flattened, Optically Thick Disk 



Ks-[24] 



FlG. 5. — Loci of colors for an optically-thick flattened disk around a MO star (black solid line) and an optically-thick "settled" disk 
around a MO star (grey solid line) as a function of disk inclination. Disks viewed edge on have extremely red colors and lie off the boundary 
of the plot. The right-most point in each plot corresponds to the face-on case. Note that since the model predicts some K-band excess, the 
colors compared to that of the stellar photosphere are redder by an additional ~ 0.1-0.2 mags. 
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Fig , 6. — SED analysis for CrA-159, comparing its emission to the razor-thin fiat optically thick disk model (e.g. Kcnvon and Ha rtmarml 
1987, dash-thrcc dots/diamonds) and a flat, optically-thick disk model constructed from the Whitney radiative transfer code (dotted line). 
For reference we overplot the stellar photosphere as a solid line. CrA-159 has disk emission consistently lying above the Whitney flat disk 
model and thus we classify it as having a primordial disk. 
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Fig. 7. — SEDs G-l (top), G-14 (middle), and G-87 (bottom). Wc fit the optical/near-IR portion of the SED to a stellar SED constructed 
from the Currie ct al. (2011) intrinsic colors (solid line), overplot the flat, optica lly-thick disk model cons t ructe d from the Whitney radiative 
transfer code (dotted line) and overplot the flat, reprocessing disk SED from Kcnyoii and Hartmann (1987, dash-three dots). All three 
stars have homologously depleted transitional disks. 
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Fig. 8. — Spectral energy distributions and histograms of inner disk radii for CrA-205, CrA-466, and G-65. For the righthand panels, the 
hatched region shows the distribution of disk inner radii from best-fit models; the grey shaded region sho ws the normalized dist ribution of 
inner radii from all models. CrA-466 and G-65 do not show evidenc e for an inner hole, in agreement with Ercolano ct al. (2009). CrA-205 
has an inner hole and thus is a transitional disk, in agreement with Sicilia-Aguilar ct al. (2008). 
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Fig. 9. — Signal-to-noise maps constructed from our reduced NGC 2362 MIPS mosaic image with the positions of ID-41 (top) and ID-63 
(bottom) circled. The color str etch in each panel go es from 1 to 4.5 a (black to blue to white). The signal-to-noise for the integrated flux 
is ~ 5 for both ID-41 and 63. Luhman et al. (2010) claim "from inspection of the mosaic" that both of these sources have SNR < 3. If 
that were the case, none of the pixels center ed on these targets would be light blue or white. However, their removal of two stars (IDs 663 
and 1091 but not IDs 809 or 931) from the[E rwin et all sample not considered here is better justified: one does has a low signal-to-noise 
and the other is blended with a different source. 
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Fig. 10. — SEDs for four homologously depleted transitional disks in Taurus, showing optical to mid-IR photometric data compared to 
the flat, optically-thick reprocessing disk limit and the lower-quartile Taurus SED. 
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Fig. 11. — SED modeling of sources shown in previous figure, displaying the best-fit SEDs from the Robitaille grid along with the 
corresponding masses for the best-fit models. For the righthand panels, the hatched region shows the normalized distribution of masses 
from best-fit models; the shaded grey region shows the normalized distribution of masses from all available models. The x-axis of the 
righthand side figures lists the disk mass, not the fractional disk mass. The division corresponding to M disfe /M* = 10' — 3) for VY Tau, 
JH-223, ZZ Tau, and FP Tau occurs at M disk = 7.5xl(T 4 , 5.5 X 1CT 4 , 4 X 10~ 4 , and 2.5 X 1(T 4 , respectively 
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Fig . 12. — K B -[8] vs. K s -[2 4l colors for K5-M2 (left panel) and M3-M6 (right panel) Taurus members with data from lLuhman et all 
(2010) and Rcbull ct al. (2010). Grey circles surrounded by squares represent transitional disks identified from SED modeling, black dots 
represent primordial disks identified from SED modeling, and black dots surrounded by squares identify disks that could either be primordial 
or transitional. SED modeling was not performed for members represented as grey dots only. Sources between the horizontal dashed line 
(K s -[8]=1.25) and the horizontal dotted line (K«- 8] 1.75) include both primordial disks and transitional disks. Those above the horizontal 
dotted line include only primordial disks while primordial disks are absent below the horizontal dashed line. 
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Fig. 13. — K s -[8] vs. fractional disk mass for homologously depleted transitional disks (grey circles surrounded by squares), "borderline" 
cases (black dots surrounded by squares), and primordial disks (black circles). The error bars identify the the interquartile range of 
disk masses for each star. For reference, we show the thr eshold in K 3 -[8 l color below which the Taurus transitional disk population begins 
(horizonatal grey line, K s -[8] = 1.75) and that assumed in Luhman ct al. (2010) for K5— M2 stars (horizontal black dotted line, K s -[8]=1.25). 
This plot shows that the ranges of disk masses for the homologously depleted transitional disk and primordial disk populations are distinct. 
The interquartile range of disk masses for the borderline cases slightly overlaps with the range for primordial disks. 



APPENDIX 

IDENTIFYING TRANSITIONAL DISKS 

Wavelength Range Required to Meaningfully Constrain Disk Masses 

To determine the wavelength range necessary to meaningfully constrain disk masses, we follow Robita ille etUl (|2007f ) 
and display SED modeling results for AA Tau, a " typical" optically-thick primordial disk (Figure [20)) . For modeling, 
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Fig. 14. — Representative SEDs for IC 348 sources. The provisional disk states for these members are (clockwise from top-left): primordial 
disk, primordial disk, homologously depleted transitional disk, and transitional disk with an inner hole. 




Fig. 15. — Representative SEDs for NGC 2362 sources. The provisional disk states for these members are (clockwise from top-left): 
primordial disk, transitional disk with inner hole, homologously depleted transitional disk, and homologously depleted transitional disk. 
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Fig. 16. — SEDs for two r\ Cha members identified as having homologously depleted transitional disks based in their photometric and 
spectroscopic flux densities (dots and thick grey line, respectively) compared to the appropriate flat, reprocessing disk models. 



we ass ume an extinction range of Ky = 1.45 to 2.45 - the average between the iFurlan et al.l (|2006l) and lRebull et al.l 
(12010D estimates ± 0.5 mags - and use the optical, 2MASS, Spitzer-IRAC/MIPS, and submillimeter data as listed by 
Robitaill e et al.1 ((2007). We adopt the same photometric errors as well. 

The best-fit SED including all these data (top-left panel) reproduces the iRobitaille et al.l results. Our distribution 
of masses (top-right) is sharply peaked about Mdisk = 0.015 M© with a median value of ~ 0.02 M©, perfectly 
overlapping with their results . Now we remove longer wavelength data to see if the predicted range of disk masses 
changes. Restricting ourselves to data equal to or shortward of 70 /jm, the dispersion in disk masses is larger, but 
clearly and sharply peaked about Mdisfc ~ 0.015 M© (bottom-left) with a median value of ~ 0.01 M©. Modeling 
optical through MIPS-24 data yields a peak noticeably shifted towards slightly lower disk masses (~ 8xl0~ 3 M©). 
More generally, the distribution of disk masses traces the shape of the distribution for all possible models (light grey 
shaded region). 

While submillimeter data is required to yield the most precise disk mass estimates, this exercise indicates that far-IR 
wavelengths effectively probe the bulk disk mass. Disk mass estimates derived from modeling optical through far-IR 
data (e.g. 70 fim) show good agreement with modeling that includes submillimeter data. Therefore, we conclude that 
a submillimeter detection is preferred but not required to estimate disk masses with a precision useful for our study. 
On the other hand, if far-IR/submm data is not available, disk mass estimates are very imprecise and likely not useful 
for our study. 

Masses of Primordial Disks and Transitional Disks 

SED modeling of optical-to-submm data for stars in 1-2 Myr-old clusters like Taurus provides useful benchmark 
values for disk properties that can be utilized to assess disk evolution. For instance, comparing near-to-mid IR emission 
from disks to the median and upper/lower-quartile Taur us SED is often used to identify sources showing evidence fo r 
warm dust depletion and/or settling in the inner disk fcf. lHartmann et alj|2005l IFurlan et aTll2006HCurrie et alJ feOOQ). 
Similarly, the median and upper/lower quartile of masses for optically-thick primordial disks compared to stellar masses 
as inferred from submillimeter measurements ("fractional disk masses") can be used to asse s s dust depletion in the 
outer disk, which traces the total mass of dust (|Beckwith and SargentHl99lt IWood et al.ll2002t lAndrews and Williams! 
2005) . Submillimeter data is particularly important in this paper since identifying homologously depleted transitional 
disks - those with low disk masses compared to primordial disks- presupposes that we can id entify appropriate lower 
limits for primordial disk masses. In this section, we revisit and expand upon the study of lAndrews and Williams] 
(|2005t ) to identify the range of fractional disk masses for primordial disks, showing that the appropriate lower limit is 
M dlsk = 0.001-0.003 M„. 

The [An drews and Willi ams! (|2005h study includes data for 153 Taurus members (their Table 1). 
lAndrews and W illiams (|2005j ) estimated disk masses for 44 of the members from SED modeling (their Table 2), 
since they have submillimeter detections and high-quality near-to-far IR photometry: we hereafter refer to this sam- 
ple as the 'Table 2 sample'. Disk masses for the other 109 members were estimated solely from their submillimeter 
fluxes/upper limits: we hereafter refer to this sample as the 'Table 1-only sample'. 

To determine fractional disk masses for disks in both samples, we adopt spectral types from iLuhman et "all (|2010T) . 
use our stellar mass estimates, and restrict the samples to K5-M6 stars. Because our goal is identifying the masses 
of optically-thick primordial disks, we further trim the sample of Class I protostars, Class III diskless stars, known 
transitional disks with inner holes (e.g. GM Aur, UX Tau), and other disks with optically-thin near-to-mid IR 
emission listed in Table [5] (e.g. V8 36 Tau, JH 223). Fin ally, we remove foregro und nonmembers (St 34) and other 
stars lacking MIPS data reported in lLuhman et al.l (|2010t ) and lRebull et al.1 (|2010h that is necessary for SED modeling 
(e.g. CIDA-11, CIDA-12, and FV Tau/c). 

Our selection criteria le aves 22 stars in the Table 1-on ly sample and 26 in the Table 2 sample. The sample of stars 
that were not modeled by lAndrews and Williams] (|2005f) . and thus listed in Table 2, is dominated by Class III diskless 
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Fig. 17. — Frequency of transitional disks vs. time. (Top panel) We compare our frequencies to the theoretical models of 
Alexander and Armitagc (2009) who determine the frequency of transitional disks vs. time from a disk evolution/planet formation model 
with a disk clearing timescale of 0.5 Myr. The dashed line corresponds to t he reference model, while the so lid line includes a dispersion 
in disk properties (initial angular momentum, disk lifetime, di s k mas s). See [Alexander and Armitagc (2009) for details. (Bottom panel) 
Our frequencies compared to those derived by[Muzcrollc ct al. (2010). Their loci correspond to frequencies assuming that all transitional 
disks have inner holes and optically thick outer regions (lower dashed line/diamonds) and allowing for a broader definition of transitional 
disk that would include many identified from this work (upper dashed line/triangles). In t his figure and the next one, the dots identify the 
average of the lower and upper limits for the transitional disk frequency as listed in Table 1121 



members. This characteristic explains the precipitous drop in the number of Table 1-only stars analyzed (22) compared 
to the nominal sample (109). 

Figure fH] shows the cumulative distribution functions for both samples, adopting lAndrews and Williams! ({2005) 
disk mass estimates and our stellar mass estimates. For the Table 2, the median and mean fractional disk masses are 
Mdisk ~ 0.011 M* and M disk ~ 0.039 M*. The interquartile range covers 0.003 M* to 0.033 M*. About 88% (23/26) 
have masses greater than or equal to 0.001 M*: the low-mass outliers are DF Tau, HN Tau, and V955 Tau. For the 
Table 1-only sample, about 82% (18/22) have masses greater than or equal to 0.001 M*: the low-mass outliers are 
CIDA-3, CoKu Tau/3, CZ Tau, and DP Tau. To more definitively dete rmine whether any member in either sample 
has a low-mass disk, we fit their SEDs using the iRobitaille et al.l ()2007l ) model grid and data reported from Spitzer 
(jLuhman et alj[2010t iRebull et al.ll2010t ) and ground-based studies, following the methods described in §4 of the main 
text. 

As listed in Table |8] and shown by Figure [22j inferred disk masses for DF Tau, HN Tau, a nd V955 Tau using the 
Robitaille models are substantially higher than values listed by lAndrews and Williams! ((2005). The median masses of 
the best-fit models are 5-9 xlO -3 M*, exceeding both the 0.001 M* limit and the previously quoted lower-quartile frac- 
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Fig. 18. — (Top) The frequency of protopl anetary disks vs. time as probed by Spitzer observations. For NG C 1333, IC 348, and 25 Ori 
the frequencies are taken from the literature (Gutcrmuth ct al. 2008; Lada et al. 2006; Hernandez et al. 2007b). The frequency for Taurus 
is taken from the number of "Class II" objects in Taurus reported bv lLuhman et all 1 120101) and the numb er of "Class III" obj ects showing 
evidence for weak IR excess (e. g. FW Tau). The frequency fo r NGC 2362 is taken from data presented in lCurrie et al.| (|200a ), which is in 
agreement with estimates from Dahm and Hillenbrand (2007). (Bottom) The transitional disk frequency vs. time for different clusters in 
our study. Overplotted are the predicted transitional disk frequencies from our parameteric model, assuming a typical protoplanetary disk 
lifetime of 3 Myr and varying amounts of that time spent as transitional disks (0.1-2 Myr). 



tional disk mass. Thus, we conclude that all of the optically-thick primordial disks modeled by I Andrews and Williamsl 
(2005|) have inferred masses greater than Mdisk = 0.001 M*. 

We obtain similar results for the low- mass outliers in our second sample. Figure [231 shows the SEDs and distribution 
of masses from best-fit Robitaille models for CIDA-3, CoKu Tau/3, CZ Tau, and DP Tau. CoKu Tau/3 and DP Tau 
have inferred masses greater than 0.001 M*. CZ Tau and CIDA-3 ha ve much lower masses and represent the only 
optically-thick primordial disks from the lAndrews and Williamsl (|2005| ) sample with inferred masses lower than 0.001 
M*. 

Figure [2H shows the revised cumulative distribution functions for the fractional masses of optically-thick primordial 
disks surrounding K5-M6 stars in T aurus. Except for the sources individually modeled in this section (listed in Table 
[S]), we adopt disk mass values from lAndrews and Williamsl (|2005f) for simplicity. Of these members, ~ 96% (46/48) 
have masses greater than or equal to 0.001 M*. The median and interquartile range for the Table 2 sample is 0.011 
M* and 0.006-0.033 M*. For the combined sample of 48 stars, these values are 5xl0~ 3 M* and 0.0024-0.015 M*; 
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Fig. 19. — The transitional disk frequency vs. time for different clusters separat ed b y spectral type: K5— M2 stars (lower locus) and 
later stars (upper locus). The error bars simply define the range in values from Table [121 taking into account disks with uncertain states. 



however, most of the masses are derived from submillimeter fluxes alone, so individual mass estimates accounting for 
the entire disk should be larger and the median and quartile disk mas ses should likewise be larger. 

Since nearly all primordial disks in the I Andrews and W illiams (2005) sample have masses greater than 0.001 M* and 
the lower quartile disk mass for the modeled sample is > 0.003 M star , we conclude that disks with masses lying below 
0.001-0.003 M* lie below the range expected for primordial disks. Thus, based on modeling current optical through 
submillimeter data, defining the lower mass limit for optically-thick primordial disks at 0.001-0.003 M* is justified. 
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Fig. 20. — Analysis of AA Tau following TRobitaille et al.l II2007D . The top-left panel shows the best-fit SED models to the AA Tau data 
from the optical through submillimeter; the top-right panel shows the histogram of disk masses from these best-fit models. The bottom two 
panels show the histogram of disk masses if we restricted our data to 70 fim or less (bottom-left panel) and 24 fim or less (bottom-right 
panel) . 
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Fig . 21. — Cu mulative distribution fu nction of primordial disk masses for 48 K5-M6 Taurus members listed in Andrews and W illiamsl 
(2005) based on Andrews and Willia ms's disk mass estimates. Also shown are distribution functions for the 'Table 1-only' sample and 
'Table 2' sample from Andrews and Williams (2005, see Appendix for definitions) identified by dotted and dashed distribution functions, 
respectively. For reference, our adopted limits for primordial disk masses - 0.001 M* and 0.003 M* — are shown as vertical solid and dotted 
grey lines, respectively. 
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Fig. 22. — SEP modeling re sults for Class II objects with optically-thic k near-to-mid IR emissi o n (pr imordial disks) in Taurus studied in 
lAndrews and Williams! J2005D that have M dlak < 0.001 M* based on the lAndrews and Williams! <|2005T ) SED modeling: DF Tau, HN Tau, 
and V955 Tau. The division corresponding to M disfc /M* = 10< - 3) for DF Tau, HN Tau, V955 Tau occurs at M diak = lxl0~ 3 , 8.5 X 
10~ 4 , and 5.5 X 10 — 4 , respectively. Using the Robitaille models, we find that all of these sources have disk masses larger than 0.001 M*. 
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Fig. 23. — SED fits and distribution of masses from best- fitting models for other primordial disks in Taurus with data from 
lAndrews and Williams! (I2005T) . The division corresponding to M disfe /M* = 10< - 3) for CIDA 3/V410 X-Ray 1, CoKu Tau/3, CZ Tau, 
and DP Tau occurs at M disk = 6 xlO -4 , 6 X 10 -4 , 5.5 X 10~ 4 , and 4 X 10~ 4 , respectively. Based on the lAndrews and Williams! 1 120051 ) 
submillimeter estimates, all of these have disk masses less than 0.001 M*. Our analysis finds that only two of the four have such low masses. 
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Fig . 24. — Cumulative distribution function of primordial disk masses for 48 K5-M6 Taurus members listed in Andrews and W illiamsl 
(2005) taking into account revised mass estimates for Taurus members discussed in the Appendix. Our results show that the disk mass 
range of 0.001-0.003 M 4 defines a meaningful lower limit for the masses of optically-thick primordial disks. 
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